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Abstract: We perform a new and updated analysis of sneutrinos as dark matter can- 
didates, in different classes of supersymmetric models. We extend previous analyses by 
studying sneutrino phenomenology for full variations of the supersymmetric parameters 
which define the various models. We first revisit the standard Minimal Supersymmet- 
ric Standard Model, concluding that sneutrinos are marginally compatible with existing 
experimental bounds, including direct detection, provided they compose a subdominant 
component of dark matter. We then study supersymmetric models with the inclusion of 
right-handed fields and lepton-number violating terms. Simple versions of the lepton- 
number-violating models do not lead to phenomenology different from the standard case 
when the neutrino mass bounds are properly included. On the contrary, models with 
right-handed fields are perfectly viable: they predict sneutrinos which are compatible with 
the current direct detection sensitivities, both as subdominant and dominant dark matter 
components. We also study the indirect detection signals for such successful models: pre- 
dictions for antiproton, antideuteron and gamma-ray fluxes are provided and compared 
with existing and future experimental sensitivities. The neutrino flux from the center of 
the Earth is also analyzed. 
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1. Introduction 

Sneutrino as a particle candidate to explain the Cold Dark Matter (CDM) present in the 
Universe has been investigated in the past in a number of interesting papers, where its relic 
abundance and its scattering cross-section off nucleons, relevant for the direct detection 
searches of dark matter, have been calculated and discussed. Sneutrinos in the Minimal 
Supersymmetric Standard Model (MSSM) have been studied in the past [Q, |3|, ^, ||, 
and from the direct detection searches they have been excluded as the major component of 
dark matter. The case of the inclusion of right-handed sneutrino fields [0, |8|, ^ 11] and 



of lepton-number violating terms [0, ^, |TT|, |l^, 13, 14, 15| have also been studied, in 

connection also to the problem of neutrino masses. It has been shown that right-handed 
components may alter significantly neutrino phenomenology, due to a reduction of the 
coupling with the Z boson. Lepton-number violating terms, which can induce a mass 
splitting of the two mass eigenstates, also lead to a modified phenomenology since they 
alter the coupling to the Z boson, which is non-diagonal in the CP basis [^] . These models 
offer a nice realization of inelastic dark matter, which has been introduced in relation to 
direct detection searches |ll|, Sneutrinos in connection with the dark matter 

problem have also been discussed in different frameworks some of which may be found in 

Refs. [H, n, n, io|, m. Hi. 



Typically, the relic abundance and the direct detection rate (which offers a very strin- 
gent experimental bound to sneutrino dark matter) have been studied and some of the 
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models have been proposed in order to circumvent the character of exclusion of sneutrino 
dark matter in the minimal version of the MSSM. However, a thorough analysis as is typ- 
ically done for neutralino dark matter, with a global study in the full parameter space of 
the supersymmetric models, has not been performed for sneutrinos. Indirect detection sig- 
nals, especially those coming from dark matter annihilation in the Galaxy, have not been 
typically discussed in the literature. 

In the present paper we wish to reconsider in a consistent way sneutrino as a cold relic 
from the early Universe and study its phenomenology relevant both for Cosmology and for 
relic-particle detection. First of all, we explicitely consider both cosmologically dominant 
and sub-dominant sneutrino configurations: in fact, we are interested not only in those 
configurations which are able to solve the CDM problem, but also those which provide a 
smaller amount of cosmological relic abundance but which could be potentially detectable 
by means of various astrophysical signals. We then study the sneutrino detection rates, both 
of direct and indirect type. Therefore, in addition to the quite relevant direct-detection 
signal, coded into the scattering cross-section off nuclei and which typically provides the 
most stringent limits to sneutrino CDM once the relic abundance bound is imposed, we 
calculate also the indirect detection signals which come from sneutrino pair-annihilation in 
the galactic environment: predictions for antiproton, antideuteron and gamma-rays signals 
are provided and discussed. This is a novel analysis for sneutrino CDM, and we will show 
that relic sneutrinos could have a chance of being detected by means of indirect searches. 
At the same time, for some particle physics models, indirect detection nicely complements 
direct searches in posing limits to the supersymmetric parameter space. Predictions for 
the neutrino flux from the center of the Earth, relevant for neutrino telescope searches, are 
also discussed. 

The particle physics models we explicitely analyze in the paper are extensions of 
the Minimal Supersymmetric Standard Model (MSSM), where the sneutrino is the scalar 
super-partner of the left-handed neutrino. Since neutrinos have masses, as is now clearly 
understood by a host of independent and very robust experimental results and theoretical 



analyses p3|, p3, 25, Eq, 27|, we will focus our attention on extensions of the MSSM which 



contain terms in the supersymmetric lagrangian which can drive neutrino masses. Connec- 
tions between neutrino physics and the phenomenology of sneutrino CDM will therefore 
arise, and we will explicitely consider them whenever relevant. The models which we will 
be considering are therefore natural and direct extensions of the MSSM which incorporate 
at the same time the new physics required to explain two basic problems of astro-particle 
physics: the origin of neutrino masses and the nature of dark matter. We do not attempt 
to be totally exhaustive on the type of supersymmetric models (something which would re- 
quire an exceedingly large analysis). Instead we concentrate on a number of the most direct 
extensions of the MSSM and derive the phenomenology of sneutrino CDM thoroughly. 

We first review the case of the Minimal Supersymmetric extension of the Standard 
Model (MSSM), by enlarging the previous analyses to the case of sub-dominant sneutrino 
CDM. This is discussed in Sect. |^. This model, not very appealing for sneutrino CDM 
and actually already almost excluded by direct detection searches (we will explicitely show 
under what circumstances a sneutrino CDM may still be viable in these model) sets the 
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basis for the following extensions. For the sake of clarity, and due to the complexity of 
the overall analysis of this paper, we will name the different classes of models with specific 
labels. The standard MSSM will be called "STD model" in the discussion. 

The first extension allows for the inclusion of a right handed field for the sneutrino 
1^, |8|, ^, 0, |lT| and is discussed in Sect. ^. We will label this class of models "LR models". 
Since we are enlarging the parameter space of the supersymmetric models by adding terms 
(and therefore parameters) to the fundamental lagrangian, we will organize our discussion 
in two steps. We first discuss the sneutrino CDM phenomenology in terms of the new added 
parameters (by keeping a fixed configuration for the rest of the parameter space, whenever 
relevant). This will be useful to disentangle specific features of the relic abundance and 
direct detection rate for this specific class of models. We then extend our analysis to 
include a full scan of the supersymmetric parameter space. All the relevant experimental 
bounds on searches of Supersymmetry at accelerators and supersymmetric contributions 
to rare processes are properly taken into account. Indirect detection signals are discussed 
for the full-scan case only. We will show that "LR models" contain viable sneutrino CDM 
candidates for a large sector of its parameter space and possess a rich sneutrino CDM 
phenomenology with potentially detectable direct and indirect detection rates. 

We then discuss a class of models with a lepton-number violating mass-term added to 



the lagrangian [12, 14, 15 1. In this case the two mass eigenstates are splitted and possess 



an off-diagonal coupling to the Z-boson which may sizeably alter both the relic abundance 



and the direct detection rate |12]. This class of models is named here "/j models" and 



is discussed in Sect. ^. In this case, the new term in the lagrangian can lead to 1-loop 
contributions to the neutrino mass, which may be potentially large depending on the values 
of the parameters involved. We therefore bound our analysis to the neutrino mass limits, 
as an additional and relevant phenomenological constraint. In this class of models, once 
all the experimental bounds are taken into account, including neutrino mass limits, we do 
not find a very interesting phenomenology for CDM: we therefore discuss only the general 
analysis, but we skip the full parameter-space scan and the discussion of indirect detection 
signals. 

The last class of models we consider contains right-handed sneutrino fields and lepton- 
flavour violating terms simultaneously 0, IC, 11, 12, 13, 14, 15] and is discussed in 
Sect. ||. These models, which we label as "MAJ models", are very interesting from the 
theoretical point of view since they can successfully accommodate both Dirac and Majorana 
neutrino mass-terms with a renormalizable supersymmetric lagrangian. Neutrino masses 
are then obtained by means of the see-saw mechanism |2|, |o|, ||, |3|, ||, |4|. A nice 
exhaustive analysis of this topic has been recently performed in Ref Also for this 

class of models we first perform a study of the parameter space relevant for the sneutrino 
sector (and by taking into account also neutrino-mass theoretical and phenomenological 
consequences). We then perform a full scan of the supersymmetric parameter space for 
two specific relevant sectors of the model, named "MAJ[A]" and "MAJ[B]", for which we 
discuss explicitely also the indirect detection rates. We will show that the class of models 
"MAJ[A]" (characterized by a TeV-scale Majorana mass-parameter in the neutrino sector) 
have again a rich phenomenology from the point of view of sneutrino CDM. 
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2. Standard minimal MSSM 



In the minimal MSSM, sneutrinos are the scalar partners of the left-handed neutrinos. 
Superfields U contains the fermionic SU(2)i doublets = (where / = e,fi,T 

runs over the three families) and its corresponding scalar doublets U = /£). The 



super symmetric lagrangian is constructed in the usual way (see, for instance, Refs. [36, p?|] ). 
We call this model "STD model". We do not explicitely detail here all the terms of the 
supersymmetric lagrangian which involve the sneutrino fields. We instead quote only the 
terms which are relevant for our discussion. The complete set of interaction lagrangians 



can be found in Ref. pq, 37, 38 1. 



The part of the superpotential relevant for the leptonic sector is (we use notations as 



in Refs. |3|] and |3|): 

W = e,,{fiH}H^ - Y/'HmR') (2.1) 

where and are the two higgs-doublet superfields, R'^ denotes the right-handed 
charged lepton superfields (which contain the right-handed sleptons) , ^ is the usual higgs- 
mixing parameters and Y^^'^ is a matrix which contains the Yukawa couplings. Repeated 
indices imply a sum-convention over them. In our analysis, which relies on a minimal 
version of the MSSM with a minimal set of relevant parameters, Y^^'^ is real and diagonal 
in flavour space and the Yukawa couplings are linked to the charged-lepton masses by 
the usual relation m/ = viY^^^ , where vi is the vacuum expectation value of the neutral 
component of the Higgs field. 

The soft-supersymmetry-breaking potential relevant for the sneutrino sector is: 

^soft = (Ml)" Lf*L/ + [e,,{Al'HmR') + h.c] (2.2) 

where A^'^ is a matrix, which we take real and diagonal in flavor space, as we do for the 
Yukawa matrix l^/"'. Also the mass matrix in the minimal version of MSSM, is taken 
to be diagonal and, in order to reduce the MSSM to a minimal set of parameters, all 
the three entries equal to a common value which we denote as Tiij^. The scalar potential 
is determined as the sum V = Vp + Vd + Kioft; where the D-term Vd describes gauge- 



interactions 37 1 and the F-term Vp is obtained by the superpotential by means of its 
derivatives over the scalar components (j)a as Vp = \dW / d(j)a\'^ ■ 

From all the above, the mass-term for the sneutrino field (for each family) is simply 
derived to be: 

^ ^ ' ^ %2Pi (2.3) 



m| + 2^-^ COS 2/3 



where P is defined as usual from the relation tan/3 = V2/V1 where V2 is the vacuum ex- 
pectation value of the neutral component of the Higgs field and mz is the Z-boson 
mass. In this case, the three sneutrinos (one for each family) are also (degenerate) mass- 
eigenstates with squared-mass m\ = m\ + 0.5 cos(2/3). Here and thereafter we will 
denote by mi the mass of the lightest sneutrino mass-eigenstate, which in order to be a 
CDM candidate must also be the Lightest Supersymmetric Particle (LSP). 

The experimental bounds on MSSM sneutrinos come from searches for supersymme- 
try at accelerators and, for sneutrinos lighter than mz/2, from their contribution to the 
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invisible Z-width: 

3/2 



Arz = — 

2 



mz J 



e{mz-2mi) (2.4) 



where Fj^ = 167 MeV is the Z-boson invisible width into one neutrino species. The bound 
we adopt is: ATz < 2 MeV |3£]. This bound constrains MSSM sneutrinos to be heavier 



than about 43.7 GeV for one sneutrino, and 44.7 for 3 degenerate sneutrino [^. but will 
be evaded in non-minimal models, as we will discuss below. 

The bound that comes from accelerator physics is induced by the non-observation 
of the corresponding charged sleptons. Contrary to MSSM sneutrinos, which are purely 
left-handed fields, charged sleptons possess both left- and right-handed components. The 
mass bound on the charged slepton depends on the assumptions made on the balance 
between left- and right-handed components (and on some more assumptions on the MSSM 
parameter space and mass splitting with neutralinos) and it is usually more conservative 
for the right-handed fields. For selectrons current limits are: 73 GeV for cr and 107 GeV 
for |4^- For smuons: 94 GeV for /l/j. For staus: 81.9 GeV for a generic mixing of 

TL and TR. 

These limits refer to the mass eigenstates of charged sleptons, which depend on two 
mass parameters: rriL, the soft-mass for the left-handed SU(2) doublet L and the soft -mass 
parameter for the right-handed sleptons m^. These parameters are in general matrices in 
flavour space, but as we discussed above we assume them as diagonal and common over 
the three families. The mass matrix for sleptons is, in this case: 

_y\^2 ^ I "T-i + "i|cos2/3(sin2 6'vi/ - ^) + (^^ + /itan/3) \ 

' I mi {Ai + fi tan j3) rnj^ — m| cos 2/3 sin^ 9w + mf I 

where m/ is the mass of the partner lepton (which is negligible, except when / = r) and 
Ai is a trilinear coupling in the soft -breaking potential. The parameter m/, is common 
with sneutrinos, while m/j is an independent parameter. Therefore, in general, a bound 
to the mass of the lightest charged sleptons does not directly transfer to a mass limit to 
the corresponding sneutrino, but depends on the relative values of and (and, to a 
lesser extent, on Ai and ^ for the r case). For instance, in a minimal Supergravity scenario 
(mSUGRA), the values of the parameters at the electroweak-scale are induced by their 
renormalization-group equation evolution from the GUT scale and they read: 

rriR = ml + 0.15 ml/2 (2-6) 
ml = mg + 0.52m^/2 (2-7) 

where mo and mi/2 are deflned at the GUT scale and are the common value for the soft 
susy-breaking mass parameters and the common gaugino mass parameter, respectively. In 
this framework the lower bound on the sneutrino mass is 84 GeV 

The version of the MSSM we adopt in this paper is a low-energy supersymmetric 
extension of the Standard Model, which does not (necessarily) invoke mSUGRA relations. 
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nij (GeV) 

Figure 1: STD model - Sncutrino relic abundance flh^ as a function of the sneutrino mass mi. The 
higgs masses for the lightest CP-even higgs is 120 GeV, for the heaviest CP-even H and the CP- 
odd A is 400 GeV. The solid (dashed) curves refer to models with (without) gaugino universality. 
The vertical line denotes the lower bound on the sneutrino mass coming from the invisible Z-width. 
The horizontal solid and dotted lines delimit the WMAP interval for cold dark matter. 

This model has been widely used to study neutralino dark matter and we extend it here to 
develop the sector where sneutrino is the LSP and the dark matter candidate. The basic 
features of the model and the parameter space not directly related to the sneutrino sector, 
together with the experimental constraints which we adopt, are reported in Appendix |^. 
In this minimal version of the MSSM, we assume that all the soft-mass parameters of the 
charged sleptons are common at the electroweak-scale, and therefore we set = rriR. In 
this case, the mass bound on the charged sleptons quoted above is transferred to a lower 
limit on the mass of the three degenerate sneutrinos which can be as low as the bound 
coming from the invisible Z width, depending on the value of tan /?. This is therefore the 
lower bound on the sneutrino mass in this class of models. 

Let us turn now to the discussion of the sneutrino phenomenology relevant for dark 
matter. First of all, we have calculated the sneutrino relic abundance, by taking into 
account all the relevant annihilation channels and co-annihilation processes which may 
arise when the sleptons are close in mass to the sneutrinos. Table |^ lists all possible 
channels for annihilation and coannihilation, referred to the first family for definiteness. 
In this minimal MSSM models, the three neutrinos are degenerate in mass: they therefore 
must be considered jointly in the calculation of the relevant processes. 

For general values of the model parameter space, sneutrinos may occasionally also co- 
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Table 1: Summary of the sneutrino annihilation and coannihilation channels. For definiteness, we 
report here the case of the first family. 
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Figure 2: STD model - Sneutrino-nucleon scattering cross section ^'''iuTicon ^® ^ function of the 
sneutrino mass mi, for the same set of parameters of Fig. |l|. The solid (dashed) curves refer 
to models with (without) gaugino universality. The vertical line denotes the lower bound on the 
sneutrino mass coming from the invisible Z-width. The dashed-dotted curve shows the DAMA/Nal 
region, compatible with the annual modulation effect observed by the experiment |45j . 



annihilate with neutralinos and/or charginos. In the models we present in this paper we 
have explicitely neglected this case, by considering configurations which possess neutralinos 
at least 30% heaver than the lightest sneutrino. Co-annihilation with neutralinos and 
charginos is more accidental than the one with sleptons, whose masses depend on some 
parameters (mostly uil) which are common with the sneutrino mass sector. The relic 
abundance is calculated by numerically solving the relevant Boltzmann equation, and the 
thermal-average of the (co) annihilation cross section is performed numerically as detailed 
in Ref. |^. We have developed a detailed numerical code for this pourpose 

An example of sneutrino relic abundance Qh'^ for the minimal MSSM is plotted in 
Fig. m as a function of the sneutrino mass (we report also the result for values of the mass 
which are below the experimental bound discussed above for the sake of the discussion in 
the following Sections). In this plot we have fixed the value of the higgs masses at 120 
GeV for the lightest CP-even state h and at 400 GeV for the heaviest CP-even H and for 
the CP-odd state A. The lightest neutralino mass is = min(294 GeV, 1.3 mi) for the 
solid curve, and = 1.3 mi for the dashed curve. This second case, which possesses light 



neutralinos, requires gaugino-non universality [^, 51, q^] in order to evade the neutralino 
mass lower bound of about 50 GeV which is instead derived for gaugino universal models. 
This is nevertheless a perfectly viable possibility. 
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Figure 3: STD model - Sneutrino-nucleon scattering cross section ^a. 



(scalar) 
clcon 



as a function of the 



sneutrino mass mi. Notations are as in Fig. H, except for the experimental curves which refer here 
to the upper limit from the CDMS experiment |4^, as re-evaluated in Ref. |Q for three 

different galactic halo models which delimit the uncertainty band. The dotted line refers to the 
standard isothermal sphere with vq — 220 km s~^ and po — 0.3 GeV cm^'^. The upper dashed line 
refers to a cored-isothermal sphere with a core radius of 5 Kpc (model Bl in Ref. [Q) and with 
Vq = 170 km s^^ and po = 0.2 GeV cm^^. The lower dashed-dotted line refers to an axisymmetric 
density profile with a power-law potential (model C3 in Ref. [Q) with vq ~ 270 km s^^ and 
Po — 1.66 GeV cm~^. 



The sneutrino relic abundance is typically very small |l|, |3|, |6|, much lower than the 



cosmological range for cold dark matter derived by the WMAP analysis |53]: 

0.092 < ftcDMh^ < 0.124 



(2.8) 



This means that sneutrinos in the minimal version of MSSM are not good dark matter 
candidates, except for masses in a narrow range which we determine to be 600-700 GeV, 
consistent with previous analyses |^]. For all the mass range from the experimental lower 
bound of about mz/2 up to 600-700 GeV sneutrinos as the LSP are cosmologically ac- 
ceptable (i.e. their relic abundance is below the WMAP bound) but they are typically 
under abundant. The three dips in the values of the relic abundance refer (from left to 
right) to the Z, h to the degenerate H and A poles in the (co) annihilation cross section, 
which occur when mi is close to half the mass of the exchanged particle. Since the mass of 
the Higgs particles are not fixed in the model, but they can span from their absolute lower 
bound of 90 GeV (which occurs for specific values of the Higgs mixing angle) to 114 GeV 
(for a SM-like light Higgs) |^, 55 1 up to several hundreds of GeV (depending on naturalness 
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assumptions), the dips shown in Fig. || may occur at different sneutrino masses. Since h 
is upper bounded in the MSSM at a value around 140-150 GeV, a dip below 75 GeV is 
always expected. The sharp drop in Qh'^ at mi ~ 80 GeV is instead due to the opening of 
the W~W^ annihilation channel. 

Dark matter direct search, which relies on the possibility to detect the recoil energy 
of a nucleus due to the elastic scattering of the dark matter particle off the nucleus of a 
low-background detector, is known to be a strong experimental constraint for sneutrino 
dark matter [^. The dependence of the direct detection rate on the DM particle rests 
into the particle mass and the scattering cross section. For sneutrinos, coherent scattering 
arises due to Z and Higgs exchange diagrams in the t-channel. The relevant cross section 
on nucleus is therefore: 

aM = <yff + o)f (2.9) 
The Z-boson exchange cross section is: 

where m_A/ denotes the nucleus mass, Af^ and Zf^ are the mass number and proton number 
of the nucleus and Gp \^ the Fermi constant. The Higgs-bosons exchange cross section is: 



47r(mi + viif^Y- 



= ..fJ^f^. + fn {A^ - ZM)f (2.11) 



where /p and /„ denote the effective coupling of the sneutrino to the nucleus, whose de- 
termination (like in the case of neutralino-nucleus scattering) is rather involved and we do 



not reproduce it here. Details may be found for instance in Refs. |56, ^]. For the case of 
sneutrinos, the effective couplings may be written as: 



i,s c,b,t 

fi = mN ( '^h + '^h 

<? Q 



I = n,p 



where kq and Kq are defined as: 



y L J 



1 



i=l,2 



27 V, ™t 

j=l,2 



with denotes the sneutrino- higgs couplings, c^, Cg are the quark- higgs couplings and /t^ 
is the nucleon mass fraction due to light quark q and /t^ = 1 — 'YTq'^'^ /t, • From the 
analyses of Ref. |5€, we derive the following values: /t„ = 0.023, /t^ = 0.034 and 



/t, = 0.14 for the proton and /t„ = 0.019, /t^ = 0.041 and /t^ = 0.14 for the neutron. We 
remind that these quantities are affected by a sizeable uncertainty which can increase the 



direct detection cross section up to a factor of a few |5C]. 
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Comparisons with experimental results are more easily and typically performed by 
using the cross section on a single nucleon ■^cr^^'jJJgQ^. Two classes of experiments are cur- 
rently running. The DAMA Collaboration uses a large-mass highly-radiopure and stable 
100-Kg Nal detector, which is specifically designed to search for the annual modulation 
effect which is expected for direct detection, as a result of the yearly relative motion of the 
Earth around the Sun 59|. The DAMA/Nal Collaboration detects an annual modula- 
tion in its low-energy rate over a period of 7 years (107731 kg x day of total exposure). 
The C.L. of this effect is 6.3 a |41, 44, When interpreted as due to a relic 

particle interaction, the allowed region shown in Fig. whose contour is outlined by the 
dot-dashed curve, arises. The allowed region shown in Fig. ^ has been calculated by taking 
into account the astrophysical uncertainties arising from galactic halo modeling [^0|. This 
is currently the only experiment which is designed to address the annual modulation effect. 
Further insight into the annual modulation effect are expected from the future results of 
the DAMA/LIBRA experiment, which is currently running with a mass of 250 Kg |61]. 

A second class of experiments do not attempt to exploit the annual modulation sig- 
nature, but instead rely on the development of background-rejection techniques in order 
to reduce the background to the sum of neutrons plus dark matter recoils. These type of 
experiments provide upper bounds on the scattering cross section. In Fig. ^ the upper 
bound from the CDMS experiment (72 kg x day with a 1.5 Kg Ge plus 0.6 Kg Si detector) 



1 46, 47 1 . For three different galactic halo models, is shown as calculated in Ref. |^8|. The 
upper experimental bound has to be regarded as a conservative limit, and the span of the 
three curves reflects the uncertainty on the galactic halo modeling. Recently results from 
the XENONIO Collaboration (15 Kg active liquid Xe) [^] have been presented (for the 
standard isothermal model only), with an exposure of 136 kg x day: these limits appear to 
be more stringent than the CDMS ones by a factor which ranges from 2 to 10, depending 
on the mass. For a recent and detailed review on direct detection experiments, see Refs. 
|6|, |4|. 

Due to the different nature of the experimental results (positive indication vs. upper 
limits) and in the absence of a solid criterion to consistently combine the various experimen- 
tal results (which are derived from different techniques and whose treatment would require 
to correctly quantify uncertainties both of theoretical origin and of experimental nature, 
which are not under our control) we will present our direct detection analysis by compar- 
ing our results separately with the DAMA/Nal region with the CDMS upper bounds (for 
which we consider the upper curve of Fig. ^ as the conservative upper limit). In order 
not to exceed with the number of figures (already very large) we will show our results by 
alternating them with either the DAMA/Nal region or the CDMS curves. The minimal 
MSSM is the only case where we present both for the same theoretical predictions. 

When compared with the DAMA/Nal annual modulation region in Fig. ^ we see that 
direct detection is indeed a strong constraint on sneutrino dark matter in the minimal 
MSSM but some very specific configurations are still viable and could explain the 
annual modulation effect. We have to consider that, whenever the dark matter particle is 
subdominant in the Universe, also its local density pQ in the Galaxy is very likely reduced 
with respect to the total dark matter density. This means that the dominant component of 
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dark matter is not the sneutrino, but still sneutrinos form a small amount of dark matter 
and may be eventually detectable. In this case we rescale the local sneutrino abundance 
by means of the usual factor: 

^ = min(l,J7/iV0-092) (2.13) 

where the value 0.092 is the lower bound for CDM abundance from WMAP, given by Eq. 
( |2.8| ). This situation is the same we encounter for neutralino dark matter, which also may 
be subdominant but nevertheless interesting to study as a relic from the early Universe. 
Since the direct detection rate is linearly dependent on the local dark matter density, and 
the experimental bounds are derived for a dominant dark matter candidate, the theoretical 
scattering cross section has to be rescaled by the factor ^. In the following, when we will 
discuss indirect detection signals, ^ will enter quadratically in the determination of the 
theoretical fluxes, since in that case the signals depend on the square of the local density. 

Comparison of Fig. |l| and ^ shows that the MSSM sneutrino is only viable when 
strongly subdominant. The mass range 600-700 GeV where the sneutrino possesses a relic 
abundance in the cosmologically relevant range is clearly excluded by direct detection. On 
the other hand, MSSM sneutrinos are still viable, but they require the special condition 
to stay on the pole of the annihilation cross section through Higgs exchange. This may 
appear a fine-tuned condition, but it is no different from the case of relic neutralinos in 
mSUGRA models, where the relic abundance is acceptable mostly in very specific regions of 
parameter space, where strong coannihilation occurs or the pole of the ^-boson exchange 
is met [ pO| , 51, 52]. We also notice that, since the masses of the Higgs bosons can vary. 



the dips in ^(^^-^^i^on ("^^^ rescaling factor) may occur at any value of the sneutrino 

mass. This means that all the mass range from the lower bound on mi up to about 200 
GeV may be compatible with the DAMA/Nal results. 



Fig. y compares our calculations with the CDMS upper limit |46, 47 1. For a standard 
isothermal halo the bound completely excludes the possibility of sneutrino dark matter. 
Some very marginal room may still be present for the conservative upper limit, when 
the sneutrino mass matches the pole condition for annihilation through a light Higgs. 
Nevertheless, this option appears to be very marginal (we remind that only the mass range 
above 70 GeV is allowed by accelerator searches). 

The results of this analysis are in agreement with the previous analyses which excluded 
the MSSM sneutrino as a dominant component of dark matter Differently from 

that studies, we explicitely consider subdominant sneutrinos as viable relic particles to be 
explored and we show that there is still an open possibility for MSSM sneutrinos, although 
very marginal. 



3. Models with a right-handed sneutrino field 

In this class of models, the neutrino/sneutrino sector is enlarged by the inclusion of right- 
handed superfields , whose scalar component are right-handed sneutrino fields 0, 
|8|, ^, 10, 11]. We call this class of models as "LR models". The right-handed fermionic 
components lead to Dirac-type masses (mj^)^ for the massive neutrinos. The relevant 
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Figure 4: LR models ~ Sneutrino mass mi as a function of the right-handed soft-mass rriN, 
for different values of the off-diagonal parameter F^: the [black] solid, [blue] dashed, [green] dot- 
dashed and [red] dot-dot-dashed refer to = 10, 10^, 10^, 10"* GeV^ respectively. The left-handed 
soft-mass is fixed here at the value of 120 GeV. 



terms in the superpotential are now: 

W = eij{^lHlH] - Yi^-^H}l]R^ + Yj-^HfUjN-^) (3.1) 

where YJ'^ is a matrix, which we again choose real and diagonal, from which the Dirac mass 
of neutrinos are obtained = V2YJ^ . Also the soft-breaking potential gets modified: 



soft 



(3.2) 

where we take both the matrices and A^"^ real and diagonal, as we do for M£ and A^"'. 
The diagonal common entries for Mfj are denoted as mj^. 

The sneutrino mass-term potential is now (we omit the family index for simplicity): 



Vrr 



m\ + ^"^i cos(2/?) + m\) 



^l^L + [m% + ml]N*N + F^{ulN + N*i'L) (3.3) 



The parameter F"^, which mixes the left- and right-handed sneutrino fields is: 

F^ = vAi, sin P — firriDCOtgP (3.4) 

In the basis defined by the vector = {u^, N*), we can define the sneutrino mass potential 
as: 



1 'J A ^2 



(3.5) 
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Figure 5: LR models ~ Sneutrino left-right mixing angle 9 a.s a function of the right-handed 
soft-mass miy, for different values of the off-diagonal parameter F^: the [black] solid, [blue] 
dashed, [green] dot-dashed and [red] dot-dot-dashed refer to = 10, 10^, 10'^, 10* GeV^, respec- 
tively. The left-handed soft-mass niL is fixed here at the value of 120 GeV. 



where the squared-mass matrix A4|^ is: 



LR 



(3.6) 



The Dirac neutrino mass is small, and can be safely neglected. The parameter mjv in 
general is independent of the other mass parameters, especially m^, which instead is related 
to the charged leptons masses, as discussed in the previous Section. We are therefore 
allowed to vary freely m^v, and whenever niN is small enough, sneutrinos lighter than 
those encountered in the previous Section are in principle viable. As long as the left- 
handed mass mi is compatible with the mass lower bounds on the charged leptons (which 
occurs for niL > 80 — 90 GeV) mi can be small, provided that mN <C uil, without entering 
in conflict with accelerator bounds. In this case, light sneutrinos may arise, and the only 
relevant limit which remains is the one provided by the invisible Z-width, which we discuss 
in a moment. 

The off-diagonal term F'^ is relevant for the mixing among the mass eigenstates, 
obtained by diagonalizing A^^,^. For A^, aligned along the Yukawa matrix Y^, i.e. for 
Ky = TjYy, F^ is necessarily very small as compared to the diagonal entries (especially the 
element {M.\^)ii of the matrix), because in that case F"^ is proportional to the neutrino 
mass rriD, and is therefore negligible since both v and ^ are electroweak-scale parameters. 
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Figure 6: LR models ~ Sneutrino left-right mixing angle 9 as a function of the right-handed 
soft-mass ttjat. Notations are as in Fig. |. The left-handed soft -mass rTii is fixed here at the value 
of 1 TeV. 



However, is in general a free parameter. In this case F"^ ~ t>A,ysin/3 may naturally be 
of the order of the other entries of the matrix, and induce a sizeable mixing of the lightest 
sneutrino in terms of left-handed and right-handed fields. We define the mixing as follows: 



Pi = — sin 9 i>L + cos ON (3 7) 

^2 = + cos 6 i>i + sin 9 N 

where 9 is the LR mixing angle. Sizeable mixings reduce the coupling to the Z-boson, 
which couples only to left-handed fields, and therefore have relevant impact on all the 
sneutrino phenomenology, as recognized in Refs. [0, ^, ^, 

The first important consequence is that the lightest sneutrino may be lighter than 
mzjl and, due to the mixing angle, still pass the invisible Z-width constraint, which now 
is modified as: 

AFz = sin^ 9 — 
2 

Also v\ annihilation and scattering cross sections which involve Z exchange are reduced 
because of the mixing. Also diagrams with higgs exchange are modified, but in a different 
way. For details. 

The free parameters in the sneutrino sector for the LR models are therefore: rrii, 
F^. We will study the model by varying mj\f and freely, while for tul will assume a 



2mi 
mz 



9{mz — 2 mi) 



(3.8) 
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Figure 7: LR models - Sneutrino mass irii vs. the sneutrino left-right mixing angle 9 for — 10^ 
GeV^ and for a scan of the soft-mass parameters tol and mjv in the ranges: 120 GeV < < 1 TeV 
and 1 GeV < m^v < 1 TeV. The rightmost full [blue] area denotes the region which is excluded by 
the invisible Z-width. The black darker area denotes the region where the sneutrino relic abundance 
is in excess of the WMAP bound; the other areas cover the region where the relic abundance is 
acceptable; in addition, in the lighter [light blue] area, the direct-detection scattering cross-section 
is inside the DAMA/Nal annual modulation region. The white areas are not covered by models for 
the values of the parameters adopted here. 



lower bound of 100 GeV, in order to assure that the mass limits on the charged sleptons is 
satisfied. 

In Figs. ^, ^ and ^ we show some features of the model relevant for our analysis. In 
Fig. ^ the mass mi of the lightest sneutrino is plotted versus the right-handed mass miy 
for various values of F'^ and for rriN = 120 GeV. Whenever rriN ^ m^, sneutrinos are 
light, and largely mixed with the right-handed component A^, as is clear from Fig. ^ where 
the sin^ is plotted against rriN, for the same set of configurations of Fig. ^. Each line 
starts from the value of ttiat below which a negative mass-squared eigenvalue occurs. 
If we increase the value of niL, the tachionic bound is met for lower values of for any 
fixed F'^, and the mixing angle diminishes (i.e. the mixing to the right-handed component 
increases). This is clear from Fig. where sin 9 vs. itln is plotted for a larger value of 
the left-handed mass parameter: mi = 1 TeV. Figs. ^ and ^ also show that, especially for 
large values of niL and F^, small mixings may be obtained also for heavy sneutrinos. 

Light sneutrinos and (very) small mixings may occur: the invisible Z-width may 
therefore be a relevant constraint. Figs. |^ and ^ show, in the plane sin^-mi and for 
different values of F"^, the region which is excluded by this constraints, as the rightmost 
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Figure 8: LR models - Sneutrino mass irii vs. the sneutrino left-right mixing angle 9 for F'^ — 10'^ 
GeV^ and for a scan of the soft-mass parameters niL and m^v in the ranges: 120 GeV < uil < 1 TeV 
and 1 GeV < < 1 TeV. Notations are as in Fig. 0, except for the lighter [light blue] area where 
now the direct-detection scattering cross-section is below the CDSM conservative upper bound. 



full [blue] area. Clearly this bound applies only for mi < mz/2, and for light sneutrinos 
it is evaded for | sm9\ < 0.4, as is clear from Eq. 3^. The scatter plot shows the distribu- 
tion of points obtained by scanning the soft -mass parameters rriL and mTv in the ranges: 
120 GeV < niL < 1 TeV and 1 GeV < mjv < 1 TeV. The off-diagonal parameter F'^ is 
fixed at different values in each Figure. We notice that a large fraction of points passes 
the Z-width bound and therefore represents viable models. We also see, in Fig. that 
for small values of F'^ the mixing angles are typically very small for light sneutrinos. Light 
eigenstates are possible only for small values of rriN, since m/, is lower bounded at 120 GeV 
in our scan, and therefore the lightest sneutrino is already mostly right-handed, which 
implies almost vanishing mixing angle. When rriN > mi the lightest sneutrino is mostly 
left-handed, but when mj\f mi, even a small F"^ can allow for rotation with small mix- 
ing angles. When F^ increases, larger mixing angles are possible also for light sneutrinos, 
and in this case the Z-width bound becomes important. This is clearly seen in Figs. ^ 
and ^ These features are quite relevant for the sneutrino relic abundance and detection 
rates, since a small mixing into vl implies reduced interactions, and therefore larger relic 
abundance and smaller direct detection rate. 



The relic abundance for the LR models is shown in Fig. 10, for a full variation of 
the sneutrino parameter space, and for a fixed configuration of the other supersymmetric 
parameters. The sneutrino parameters are varied in the ranges: 120 GeV < m^ < 1 TeV, 
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Figure 9: LR models - Sneutrino mass irii vs. the sneutrino left-right mixing angle 9 for — 10** 
GeV^ and for a scan of the soft-mass parameters niL and m^v in the ranges: 120 GeV < uil < 1 TeV 
and 1 GeV < < 1 TeV. Notations are as in Fig. 0, except for the lighter [light blue] area where 
now the direct-detection scattering cross-section is below the CDSM conservative upper bound. 

1 GeV < rriN < I TeV and 10 GeV^ < < 10^ GeV^. The other parameters are 
fixed at the same values used in the previous Section for Figs. ||, || and § the lightest 
neutralino is 30% heavier than the sneutrino (which implies gaugino non-universality when 
the neutralino mass is light) and higgs masses of 120 GeV for the lightest CP-even higgs 
h and 400 GeV for the heaviest CP-even H and the CP-odd A. All the models shown in 
the plot are acceptable from the point of view of all experimental constraints, including 
the invisible Z-width. The horizontal solid and dotted lines delimit the WMAP interval 
for cold dark matter. 

Contrary to the previous case of the minimal MSSM model, in LR models sneutrinos 
may represent the dominant dark matter component for a wide mass range, which extends 
(for the specific supersymmetric configurations discussed here) from a few GeV to about 
300 GeV. This is in fact due to the reduced Z coupling which occurs as a consequence of 
the mixing to the right-handed field A^. The dips in the scatter plot are again due to the 
Z, h and {H,A) poles in the annihilation cross section (from left to right). We remind 
that the location of the higgs poles is a consequence of the choice of the higgs masses 
we are using for this specific case. The sharp drop at mi = mw is due to the opening 
of the annihilation channel into the W~^W~ pair. The points with an acceptable relic 
abundance at masses below 10 GeV are obtained for light neutralinos, since in this case the 
annihilation cross section into neutrinos through the exchange of neutralinos gets enhanced, 
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Figure 10: LR models - Scatter plot of the sneutrino relic abundance flh^ as a function of the 
sneutrino mass mi, for supersymmetric models with the lightest neutralino 30% heavier than the 
sneutrino (which implies gaugino non-universality when the neutralino mass is light) and higgs 
masses of 120 GeV for the lightest CP-even higgs h and 400 GeV for the heaviest CP-even H 
and the CP-odd A. The sneutrino parameters are varied as follows: 120 GeV < tol < 1 TeV, 
1 GeV < ruN < 1 TeV and 10 GeV^ < < lO'' GeV^ All the models shown in the plot are 
acceptable from the point of view of all experimental constraints. The horizontal solid and dotted 
lines delimit the WMAP interval for cold dark matter. 



and therefore reduces the relic abundance to acceptable levels. This possibility occurs if 
the gaugino sectors possesses non-universality [ |6^ ] (otherwise the lightest neutralino mass 
is lower bounded at the value of 50 GeV by LEP2 analyses H, 54, |6^): gaugino non- 
universality therefore appears instrumental in allowing light sneutrinos to be cosmologically 
viable ||, ^jT^, 

Figs. 1^ ^ and ^ show how the configurations with acceptable relic abundance are 
distributed in the sin 9-mi plane. When is small, it typically induces very small mixings 
for light sneutrinos. This is clearly shown in Fig. |^ which refers to = 10^ GeV^. Light 
sneutrinos are possible in these models, and due to the small mixing with the Z boson they 
easily evade the invisible Z-width bound. However, the small-mixing angle, typically much 
smaller than 2 • 10~^, leads to values of the relic abundance in excess of the cosmological 
bound. The strong suppression of the Z coupling for sneutrinos lighter than myy in these 
models with small -F^, require the opening of the W~^W~ channel in order to reduce the 
relic abundance to acceptable levels. Only for sneutrinos in the mass range 80-200 GeV 
cosmologically acceptable configurations are found. When F"^ increases, the cosmologically 
allowed parameter space opens up. For F^ = 10^ GeV^, shown in Fig. |, mixing angles are 
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Sneutrino-nucleon scattering cross section ^cr|^^'j.^[|j^'^ as a function of the 



Figure 11: LR models 

sneutrino mass mi, for = 10^ GeV^ and for a scan of the soft-mass parameters m^, and tojv in 
the ranges: 120 GeV < < 1 TeV and 1 GeV < ttitv < 1 TeV. The other model parameters are as 
in Fig. The solid curve shows the DAMA/Nal region, compatible with the annual modulation 
effect observed by the experiment ^ All the configurations shown possess relic 

abundance compatible with the WMAP bound. 



typically larger than in the previous case also for light sneutrinos, and they may reach values 
of 0.2. Cosmologically acceptable configurations are obtained for masses which range from 
10 GeV up to 200-300 GeV. Larger values of may allow for a further enlargement of 
the allowed mass range: for = 10^ GeV^, shown in Fig. ^, sneutrinos are cosmologically 
allowed starting from a few GeV mass up to 200-300 GeV. 

The sneutrino-nucleon cross section C'^nu'cfeon' same super symmetric configu- 

rations used in Fig. 10, is shown in Figs. ^ and 13 for progressively larger values of 
F"^. The direct detection bound, although relevant for many configurations, now is easily 
evaded. Most of the configurations are allowed, and a large fraction of these would be a 
candidate to explain the DAMA/Nal annual modulation effect. The compatibility with 
sneutrino dark matter and the annual modulation effect increase when F'^ is increased, 
especially for light sneutrinos. Fig. 11, which refers to F'^ = 10^ GeV^ shows the lower 
bound of 80 GeV in the sneutrino mass which is due to the relic abundance bound, and that 
was discussed in relation to Fig. ^. For cosmologically acceptable sneutrinos, most of the 
configurations are compatible with direct detection searches. Larger values of F'^ open up 
the possibility of lighter sneutrinos, as is shown in Fig. 12, where F^ = 10^ GeV^ has been 
adopted. All configurations for sneutrinos lighter than 60 GeV are currently compatible 
with direct detection searches. Fig. 13, which refers to F'^ = GeV, i.e. to a value of 
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Figure 12: LR models - Sneutrino-nucleon scattering cross section ^'''nuTicon ^ function of the 
sneutrino mass mi, for = 10^ GeV^, 120 GeV < thl < 1 TeV and 1 GeV < mjv < 1 TeV. The 
other model parameters are as in Fig. nO. The dashed, dotted and dot-dashed curves denote the 
CDMS upper bounds ||, as in Fig^. 



the same order of magnitude of the diagonal mass terms in M\ji of Eq. show that 

in the whole mass range from few GeV up to 200 GeV sneutrinos are compatible with the 
annual modulation effect, without invoking any fine-tuned condition on the parameters. 
This is a quite remarkable feature of LR models. 

The analysis presented so far, useful to discuss the features of the sneutrino relic abun- 
dance and direct detection rate, was specific to fixed values of the higgs masses and peculiar 
values of the neutralino and chargino masses. We now extend our analysis to a full scan of 
the supersymmetric parameter space, analogous to the studies for neutralino dark matter. 
The parameter space is defined in Appendix ^, where also the experimental constraints on 
the supersymmetric model are discussed. Parameters specific to the sneutrino sector are 
varied in the following intervals: 100 GeV < m/, < 3 TeV, 1 GeV < ra^ < 1 TeV and 
1 GeV^ < f2 < 10^ GeV^. 



The sneutrino relic abundance is shown in Fig. 14. The most relevant new feature 
is that for the full supersymmetric scan, the mass range allowed by the cosmological con- 
straints is enlarged up to 800 GeV, and all the mass interval above the Z-pole may lead 
to strongly subdominant sneutrinos. This is due to either to the occurrence of the higgs 
poles in the annihilation cross section discussed above or to the mixing with the right- 
handed field. The occurrence of sizeable mixings with the N field is specially important 
for allowing lighter sneutrinos and for enhancing the relic abundance around the Z pole. 
Prom Pig. 14 we can conclude that, for a full scan of the supersymmetric parameter space 
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Figure 13: LR models - Sneutrino-nucleon scattering cross section ^'''iticloon ^ function of the 
sneutrino mass mi, for = 10'' GeV^, 120 GeV < < 1 TeV and 1 GeV < mjv < 1 TeV. The 
other model parameters are as in Fig. |l^. The solid curve shows the DAMA/Nal region, compatible 
with the annual modulation effect observed by the experiment O, li^, 44, Esl. 



in LR models, after all experimental (and theoretical) constraints are imposed, sneutrino 
dark matter is perfectly viable, both as a dominant and as a subdominant component, for 
the whole mass range 15 GeV < mi < 800 GeV. The lower limit of 15 GeV represents 
therefore a cosmological bound on the sneutrino mass in LR models, under the assumption 
of i?-parity conservation. This result may be confronted with the one obtained for relic 
neutralinos in gaugino non-universal models |]66| , [7C| , [7l| , 72|, which is about 6 GeV. 



The sneutrino-nucleon cross section is shown in Fig. |15|. Only points which are 
accepted by the cosmological constraint are shown. We see that the presence of the mixing 
with the right-handed N fields opens up the possibility to have viable sneutrino cold 
dark matter. A fraction of the configurations are excluded by direct detection, but now, 
contrary to the minimal MSSM case, a large portion of the supersymmetric parameter 
space is compatible with the direct detection bound, both for cosmologically dominant 
(denoted by [red] crosses) and subdominant ([blue] points) sneutrinos. The occurrence of 
sneutrinos which are not in conflict with direct detection limits and, at the same time, are 
the dominant dark matter component, is a very interesting feature of this class of models. 
Fig. m compares our theoretical calculations with the DAMA/Nal region (comparison with 
the CDMS upper bounds in straightforward): it is clearly seen that sneutrino dark matter 
could explain the annual modulation effect (as well as neutralinos do in many realizations 
of supersymmetric theories [73, 74, [f5|]). 
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Figure 14: LR models - Scatter plot of the sneutrino relic abundance flh^ as a function of the 
sneutrino mass mi, for a full scan of the supersymmetric parameter space, as explained in the text. 
The sneutrino parameters are varied as follows: 100 GeV < ttt-l < 3 TeV, 1 GeV < ttiat < 1 TeV 
and 1 GeV^ < F"^ < 10^ GeV^. The horizontal solid and dotted lines delimit the WMAP interval 
for cold dark matter. 



The correlation between the direct-detection relevant cross section ■^'^^^'^c^/pon^nd the 
relic abundance is shown in Fig. 16. Cosmologically dominant (or slightly subdominant) 
sneutrinos are compatible with the current level of sensitivity of direct-detection experi- 
ments, as discussed above. Nevertheless, a fraction of these configurations refer to direct- 
detection cross-sections which are up to 4 orders of magnitude below current sensitivities. 
At the same time, a fraction of the configurations which are at the level of direct detection 
sensitivity refer to quite small values of the sneutrino relic abundance: these subdominant 
sneutrinos correspond to a case where the dark matter is mostly composed of a different 
candidate, but nevertheless represent a relic from the early Universe potentially detectable 



in the laboratory (a very interesting situation by itself!). Fig. |16| shows that, provided 
a suitable model like the LR one discussed here, sneutrino dark matter may exhibit a 
phenomenology which is not very different from the neutralino in a low-energy effective 
MSSM. 
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Figure 15: LR models - Sneutrino-nucleon scattering cross section C^nucicon ^ function of the 
sneutrino mass mi, for a full scan of the supersymmetric parameter space. Parameters are varied 
as m Fig. ||. [Red] crosses refer to models with sneutrino relic abundance in the cosmologi- 
cally relevant range; [blue] open circles refer to cosmologically subdominant sneutrinos. The solid 
curve shows the DAMA/Nal region, compatible with the annual modulation effect observed by the 
experiment |l|, ||, ||, ||, ||. 



Figs. 17, 18 and show the distribution of the allowed configurations in three different 
sections of the sneutrino parameter space. In all three figures, [red] crosses refer to models 
with sneutrino relic abundance in the cosmologically relevant range; [blue] dots refer to 
cosmologically subdominant sneutrinos; open dots mark the configurations which have a 
direct-detection cross section in the current sensitivity range. Fig. 17 shows that, when mjv 
is below 100 GeV, cosmologically allowed configurations require mi, to be close to its lower 
limit: this is necessary in order to have a large enough mixing angle 9, which otherwise 
would produce an exceedingly large relic abundance. The off-diagonal parameter needs 
also to be tuned accordingly, in a range from 10^ GeV^ to 10^ GeV^, as displayed in Fig. 



18. When m^v crosses the lower bound on tul, more wide possibilities open up. Light 



sneutrinos, with masses below 40-50 GeV, need values of the mixing angle such that sin 6 
is around 0.05-0.5, with the specific correlation shown in Fig. 19. The correlation with 
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Figure 16: LR models - Sneutrino-nucleon scattering cross section ■fcr^j^^^j^^jj vs. the sneutrino relic 
abundance flh^, for a full scan of the supersymmetric parameter space. Parameters are varied as in 
Fig. The horizontal [green] band denotes the current sensitivity of direct detection experiments; 
the vertical [yellow] band delimits the WMAP interval for cold dark matter. 



values of Cc"nucieon^^ Current direct-detection sensitivity range are also displayed: the 
most relevant feature is that they require F"^ to be large when mjv is large, in order to 
procude enough mixing to reduce the direct detection cross-section to acceptable levels at 
large sneutrino masses. 

Now that we have assessed the possibility to have viable sneutrino dark matter candi- 
dates, let us now move to the study of their indirect detection signals. Dark matter particles 
distributes in the galactic halo may annihilate in pairs and produce a bunch of possible 
signals. Among these signals antimatter, namely antiprotons and antideuterons ]|7^, and 
gamma-rays may be produced. Exhaustive and detailed studies have been performed for 
the case of neutralino dark matter, but not yet for sneutrino relics. 

We start by studying the antiproton signal. The antiproton source term qp^{r, z]Tp) 
at a given position in the Galaxy, defined by the radial component along the galactic plane 



-25- 



T 1 1 1 1 I I I I, y 1 1 1 1 — I I I I I . 




Figure 17: LR models - Distribution of cosmologically acceptable models in the itin - rriL plane, 
when a full scan of the supersymmetric parameter space is performed. Parameters are varied as 
in Fig. [Red] crosses refer to models with sneutrino relic abundance in the cosmologically 

relevant range; [blue] dots refer to cosmologically subdominant sneutrinos; open dots mark the 
configurations which have a direct-detection cross section in the current sensitivity range. 



r and the vertical component z, is defined as: 



r(o-annf)o 9p{Tp) 



nil 



(3.9) 



where Tp is the antiproton kinetic energy, (cTanni')o denotes the sneutrino annihilation cross 
section averaged over the velocity distribution of sneutrinos in the Galaxy (which are 
strongly non-relativistic since their average velocity is of the order of /? ~ 10^^), gpiTp) is 
the antiproton energy-spectrum per annihilation event and pi,-^ (r, z) denotes the sneutrino 
density distribution, which is assumed to be proportional to the total dark matter density 
distribution as py^{r,z) = ^pdm('^, -z), in order to take into account both dominant and 
subdominant sneutrino relics. The energy-spectrum gp{Tp) is a sum over the different 
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Figure 18: LR models - Distribution of cosmologically acceptable models in the niN ~ plane, 
when a full scan of the supersymmetric parameter space is performed. Parameters are varied as 
in Fig. |lj. [Red] crosses refer to models with sneutrino relic abundance in the cosmologically 
relevant range; [blue] dots refer to cosmologically subdominant sneutrinos; open dots mark the 
configurations which have a direct-detection cross section in the current sensitivity range. 



energy spectra produced by the various final states of the annihilation process: 

9piTp) = J2 - F) i-^j (3.10) 

where F labels the final states. We have modelled the spectra dN? /dTp as discussed in 
Refs. [^], by means of a semi-analytic calculation which follows the production and decay 
chain of each final state until a quark is produced. We use detailed fits and interpolations 



over the results of PYTHIA [78] simulations for the treatment of the processes involved in 



the quark hadronization the subsequent hadron decays [77]. The antiproton source spectra 
are then propagated in the galactic environment to determine the antiproton fiux at the 
local position in the Galaxy <I>(i?07 0, Tp): 



q^^{r, z- Tp) ^Ro, 0, Tp) (3.11) 
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Figure 19: LR models - Distribution of cosmologically acceptable models in the mi - sin 6* plane, 
when a full scan of the supersymmetric parameter space is performed. Parameters are varied as 
in Fig. [Red] crosses refer to models with sneutrino relic abundance in the cosmologically 

relevant range; [blue] dots refer to cosmologically subdominant sneutrinos; open dots mark the 
configurations which have a direct-detection cross section in the current sensitivity range. 



We model the galactic environment as a two-zone diffusion model and for the solution 
of the propagation equation, which takes into account antiproton diffusion, scattering, 
annihilation, energy losses, propagation against the galactic wind and reacceleration, we 
use the results of Ref. [77|, where a detailed analysis which takes into account astrophysical 
uncertainties (relevant for the antiproton signal) are discussed. Theoretical uncertainties 
are large and we will comment on this point: in all our figures we will show the result for 
the median estimate of the antiproton signal, as defined in Ref. [77]. The antiproton flux 
has then to be propagated against the solar wind in order to provide fluxes at the Earth 
position, i.e. " top-of-atmosphere" (TOA) fluxes. We will show our results for a period of 
solar minimum activity. 

Antiproton fluxes from sneutrino annihilation in the galactic halo are provided in Fig. 
|20| and |2l| which show scatter plots of fluxes calculated at fixed antiproton kinetic energies: 
Tp = 0.23 GeV in Fig. |o| and Tp = 37.5 GeV in Fig. ||. The scatter plots refer to 
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cosmologically acceptable configurations: the [red] crosses denote cosmologically dominant 
sneutrinos, the [blue] points refer to the case of subdominant sneutrinos. The small grey 



points show those configurations which are excluded by direct detection. In Fig. 20 



which refers to the flux in a low-energy bin where the antiproton signal may have better 
possibility to be disentangled by the background which is due to cosmic-rays spallation over 
the galactic medium. The [yellow] shaded area denotes the amount of exotic antiprotons 
which can be accommodated in the BESS data |7£, ^ in that energy bin. This has 
been established on the basis of the theoretical calculation of the antiproton background 
1 81] and of the BESS measurements [79, |8^, by determining the maximal amount of exotic 
antiproton flux which can be accommodated on the top of the background, without entering 
in conflict with the BESS data and its experimental error, at 90% C.L. We see that the 
theoretical predictions are not currently excluded by BESS: the maximal predictions, which 
occur for sneutrino masses in the range 50-200 GeV, are at least one order of magnitude 
below the current BESS bound. However, we have to remind that the theoretical estimates 
of the antiproton signal are largely affected by astrophysical uncertainties related to the 
knowledge of the parameters which enter the diffusion equation. It has been shown in Ref. 
1 77] that this uncertainty can lead to an increase of about a factor of 8 or a decrease of 
up to a factor of 10. Therefore, our prediction in Fig. 2C may be altered by this factor 
for different choices of the propagation parameters in their allowed ranges [82]. In the case 
of the choice of astrophysical parameters which produce the maximal antiproton signal, 
the scatter plot in Fig. ^ would be enhanced by a factor of 8: in this case, still BESS 
data would not exclude any sneutrino configuration, but all the mass range from 50 to 
200 GeV would have configurations potentially detectable with just a small increase in the 
experimental sensitivity. 

Prospects for the future are shown by the dashed and dotted horizontal lines, which 
denote our estimated sensitivity of the PAMELA |8|] (dashed hue) and AMS [Q (dotted 
line) detectors to exotic antiprotons after a run of 3 years: the sensitivities are determined 
as admissible excess within the statistical experimental uncertainty if the measured antipro- 
ton flux consists only in the background (secondary) component. The estimate has been 
performed by using the background calculation of Ref. 1 81 ] , and refers to a 1-a statistical 
uncertainty. All the supersymmetric configurations in Fig. 20 above the dashed or dotted 
lines can be potentially identified by PAMELA or AMS as a signal over the secondaries, 
while those which are below the dashed or dotted lines will not contribute enough to the 
total flux in order to be disentangled from the background. We therefore see that, in the 
case of the median antiproton flux shown in the Figure, only AMS will have the possibility 
to detect a signal from sneutrino dark matter, for masses around 60 GeV or for masses in 
the range 100-200 GeV. We notice that most of these configurations refer to subdominant 
sneutrinos: this justifies the approach to consider relic particle candidates also when they 
are not dominant dark matter components, since a signal from relic particles in the galaxy 
may well be discovered. 

Fig. 20 also shows that, in the case of astrophysical propagation parameters close to 
the values which provide the maximal antiproton signal, both PAMELA and AMS will 
have good chances of detection, for sneutrinos in the mass range 50-200 GeV. We also 



-29- 



notice that in this case a large fraction of configurations with masses in the range 65-130 
GeV, and which could be potentially detectable by AMS, are actually already excluded 
by direct detection (grey points). This shows a very nice interplay between different dark 
matter searches techniques. This is further elucidated in Fig. 22, where the antiproton 
signal is plotted against the direct detection cross section (.(^^^^leol- '^^^ current bound 
from both BESS data and from direct detection experiments is shown, as well as the 
foreseen capabilities of PAMELA and AMS, together with the current direct-detection 
sensitivity region. From this figure we see that a fraction of the configurations which are 
currently excluded by direct detection would have been in reach of PAMELA and AMS. 
More interestingly, Fig. 22 shows that direct detection and antiproton searches offer a 
good deal of complementarity. Direct detection is sensitive to configurations which will 
be hardly probed by antiproton searches (those points inside the vertical [green] band and 
which refer to low antiproton fluxes). On the contrary, some configurations which refer to 
a large antiproton signal but are below current direct detection sensitivity are also present. 
Nevertheless, from this figure direct detection appears to be a stronger probe to sneutrino 
dark matter (in the sense that it can explore a wider region of the parameter space) but 
good chances of antiproton detection are also present. 

We come back now to Fig. |2l|, where the antiproton signal is calculated for a higher 
antiproton kinetic energy: Tp = 37.5 GeV. At this energy, the CAPRICE experiments 
reports the detection of a flux which is potentially in excess of the theoretical background 
||85|] . Although this excess has a small statistically significance, nevertheless is an intriguing 
possibility to be studied. In Fig. ^ we show the band which refers to a signal which would 
fill the CAPRICE excess. The scatter plot of the theoretical prediction for sneutrino 
dark matter is not able to reach the level of the CAPRICE excess, for the median choice 
of astrophysical parameters. We nevertheless comment that for the maximal choice, a 
marginal compatibility would arise, for sneutrino masses in the range 200-500 GeV. In 
Fig. ^ we also show our estimated sensitivities for PAMELA and AMS, determined as we 
discussed above for the lower energy bin. For the median choice of astrophysical parameters, 
AMS will have a marginal detection potential, while for the maximal case, both AMS and 
PAMELA will probe configurations in the mass range 200-500 GeV. Light sneutrinos are 
not probed at this antiproton energies: they cannot produce antiprotons at energies above 
their mass, since annihilation occurs almost at rest. 

Another indirect detection signal which is very promising is the production of an- 
tideuterons |76]. In this paper especially it has been shown that the low energy tail of the 
antideuteron signal offers a good possibility to disentangle the signal from the background, 
since kinematical conditions allow low-energy antideuterons to be produced in the dark 
matter pair annihilation easily, while the spallation process is suppressed below a kinetic 
energy per-nucleon of = 1-3 GeV/n. The predictions for sneutrino dark matter are 
shown in Fig. 23 for Tg = 0.23 GeV/n. No experimental limit is currently available, 
but detectors, which will be able to reach good sensitivities, are under development: Fig. 



2| shows the expected sensitivity of the GAPS @, ^ and AMS [g] detectors. We see 



that antideuteron searches will offer a sensitivity to sneutrino dark matter similar to an- 
tiproton searches. A signal detectable in one channel by two different detectors, will be 
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Figure 20: LR models - Antiproton flux at the antiproton kinetic energy Tp — 0.23 GeV as 
a function of the sneutrino mass rni, for the galactic propagation parameters which provide the 
median value of antiproton flux and for a solar activity at its minimum. [Red] crosses refer to 
models with sneutrino relic abundance in the cosmologically relevant range; [blue] dots refer to 
cosmologically subdominant sneutrinos; light gray points denote configurations which are excluded 
by direct detection searches. The [yellow] shaded area denotes the amount of exotic antiprotons 



which can be accommodated in the BESS data |79, pOl]. The dashed and dotted lines show the 



PAMELA |83l and AMS pS sensitivities to exotic antiprotons for 3 years missions, respectively. 



detectable also in the other channels, again in two different experiments. Moreover, one of 
the detector, AMS, has the capabilities to detect a signal in both channels. The possibility 
to cross-correlate different signals and to complement their information would be an ex- 
traordinary opportunity for dark matter searches. This is further complemented by direct 
detection, as discussed above. 

We now move to discuss the signal which consists in the production of a diffuse gamma- 
rays flux. Sneutrino annihilations may produce, through the decay chain of their annihila- 
tion products, also gamma-rays, which mostly come from the decay of neutral pions and 
other mesons produced in the annihilation process. The gamma ray flux arriving at the 
Earth from a given angular position ip in the sky is the integral along the line of sight of 
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Figure 21: LR models - Antiproton flux at the antiproton kinetic energy Tp — 37.5 GeV as a 
function of the sneutrino mass mi. Notations are as in Fig. ^ except for the upper [yellow] shaded 
band, which delimits the possible excess over the background in the CAPRICE data The lower 
[yellow] area refers to fluxes compatible with CAPRICE. 



all the differential contributions of sneutrino annihilation: 

^r(^7,^) = 9,{E,) m (3.12) 



where g^{E^) is the gamma-ray spectrum, defined analogously as we defined above the 
antiproton spectrum gp{Tp) and I{tp) denotes the line-of-sight integral of the squared of 
the dark matter density in the direction ip: 



IW= / />ijr(A,^)]dA (3.13) 

^l.o.s. 



For the production of gamma-rays arising from the hadronization of quarks and the sub- 



sequent decay of hadrons, we use the detailed fit of PHYTIA simulations of Ref. |9£, 97 1. 
Fig. 24 shows our predictions for the gamma-ray flux at = 1.5 GeV coming from 



the center of the Galaxy, in an angular bin which corresponds to the EGRET |9C, ^] field 
of view. The choice of the energy bin refers to the case where EGRET detects an excess 
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Figure 22: LR models - Antiproton flux at the antiproton kinetic energy Tp — 0.23 GeV vs. the 
sneutrino-nucleon scattering cross section '?o'jj^'j!^|,"j|. [Red] crosses refer to models with sneutrino 
relic abundance in the cosmologically relevant range; [blue] dots refer to cosmologically subdominant 



sneutrinos. The horizontal solid line denotes the upper limit from BESS \79, pO| and the [yellow] 



shaded area shows the amount of exotic antiprotons which can be accommodated in the BESS 
data. The dashed and dotted horizontal lines show the PAMELA [|3) and AMS [Q sensitivities to 
exotic antiprotons for 3 years missions, respectively. The vertical solid line denotes a conservative 
upper limit from direct detection searches and the [green] vertical shaded area refers to the current 
sensitivity in direct detection searches. 



of gamma-rays over the background from the galactic center. The [yellow] shaded area in 
facts refers to this excess: exotic gamma-ray fluxes inside this band are compatible with 
the EGRET measurement and those close to the solid horizontal line, which delimits the 
area, are able to explain the excess. 

The gamma-ray signal is strongly sensitive to the behavior of the dark matter density 
profile toward the inner regions of the Galaxy. In Fig. ^ we have used a strongly peaked 
profile: the radial behavior is r~^'^, like the Moore et al. profile ]^8| , 8£] or those obtained 



from milder distributions by effects due to the growth of a black hole 

-1 



baryon dissipation [101, 102]. For a NEW profile with a 



I, H, |100| or of 
behavior, our theoretical 
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Figure 23: LR models - Antideuteron flux at the antideuteron kinetic energy (per nuclcon) Tp — 
0.23 GeV/n, as a function of the sneutrino mass mi. Notations are as in Fig. The dashed and 
dotted hnes show the GAPS IstII and AMS Hi sensitivities. 



estimates would decrease by a factor of 60 and an additional factor of 10 for cored profiles 
|9|. 

Fig. 24 shows that, for an r^^-^ profile, sneutrino dark matter is at the level of explain- 
ing the EGRET excess, for masses in the range 60-250 GeV. Most of these configurations 
refer to subdominant sneutrinos. Clearly for an NFW profile the predicted fluxes are too 
low to fill the EGRET excess. 



Fig. 24 also shows our estimate for the capabilities of GLAST, for a 1-year data- 
taking. We have taken into account GLAST effective area as shown in Ref. |92|, and 
we have derived our predictions for the same angular energy bin of EGRET. We see that 
GLAST will be sensitive to configurations of masses between 30 GeV and 600 GeV, and will 
be close to access also very light sneutrinos with further live-time of data. Again we have 
to warn that in the case of less steep profiles, the sensitivity of GLAST will be restricted to 
a smaller mass range (60-300 GeV), but nevertheless a good deal of configurations will be 
probed. We comment also that the angular resolution of GLAST will be much better than 
the EGRET one and therefore the capabilities of GLAST would be even more promising 
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Figure 24: LR models - Gamma-ray flux from the galactic center at the photon energy — f .5 
GeV, as a function of the sneutrino mass mi, for a galactic profile of Moore's type ||8^, ^ and 
for the angular resolution of EGRET |90[ |9^. [Red] crosses refer to models with sneutrino relic 
abundance in the cosmologically relevant range; [blue] dots refer to cosmologically subdominant 
sneutrinos; light gray points denote configurations which are excluded by direct detection searches. 
The [yellow] shaded area denotes the amount of exotic gamma-rays compatible with the EGRET 
excess [^0[ ^ . The dashed line shows the GL AST ]Q sensitivity for a 1 year data-taking and for 
the same EGRET angular bin. 



than those shown in Fig. 

Finally, we mention that indirect detection of dark matter could also be performed at 
neutrino telescopes, since dark matter may accumulate in the central regions of the Earth 
and the Sun by gravitational capture and there annihilate. The only annihilation products 
which can escape are neutrinos, and these can be searched for as upgoing muons in a 



neutrino telescope. In Fig. 25 we show the predictions for upgoing muons from the Earth 
1 105 ] and compare them with the current experimental limits from SuperKamionande, 
MACRO and AMANDA. We see that neutrino telescopes are sensitive to a large fraction 
of sneutrino configurations, although many of the configuration which are excluded by this 



technique are also excluded by direct detection. This is manifest in Fig. 26, which shows 
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Figure 25: LR models - Upgoing muon flux from sncutrino pair anniliilation in the center of 
the Earth $Earth^ ^ function of the sncutrino mass mi. [Red] crosses refer to models with 
sneutrino relic abundance in the cosmologically relevant range; [blue] dots refer to cosmologically 
subdominant sneutrinos. The solid, dashed and dot-dashed lines denote the upper limits from the 
SuperKamiokande [||, MACRO [Q and AMANDA [|| detectors, respectively. 

the correlation of the upgoing muon flux from the Earth with the direct detection cross 

. . J. (scalar) 

section eCcieon- 

4. Models with a lepton— number violating term 

Models with lepton-number violating terms can allow for Majorana neutrino masses. 
The most direct way to include a Majorana mass term is to introduce a gauge-invariant 
dimension-5 operator of the type [^, 13, 35, 



^ = j^i'^jLiH,)ieMLiHi) + h.c. (4.1) 

In this case, a Majorana mass term for the neutrino is generated when the neutral compo- 
nent of the Higgs field acquires a vacuum expectation value and the neutrino mass which 
arises is of the order of tjim ~ gv'^ /M\. This can be made compatible with neutrino mass 
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Figure 26: LR models - Upgoing muon flux from sneutrino pair anniliilation in the center of 
the Earth (J)^^''*'' vs. the sneutrino-nucleon scattering cross section ^'''iuTicon- [R-^d] crosses refer 
to models with sneutrino relic abundance in the cosmologically relevant range; [blue] dots refer to 
cosmologically subdominant sneutrinos. The horizontal solid line denotes the current upper limit 
from neutrino telescopes. The vertical solid line denotes a conservative upper limit from direct 
detection searches and the [green] vertical shaded area refers to the current sensitivity in direct 
detection searches. 



bounds for Ma close to the GUT scale. The dimension-5 operator is clearly not funda- 
mental as an extension, and makes the supersymmetric lagrangian non-renormalizable. 
Nevertheless, this operator may arise as an effective term from new physics at the high 
energy scale M\. 

If we apply this extension at the MSSM lagrangian, we are allowed for a L-number 
violating term also in the sneutrino lagrangian, which modifies the mass-term potential as: 

'>1'>L + \m%{9Li>L + (4.2) 
where is a mass parameter that makes the mass lagrangian no longer diagonal in the 



m-L + 2"^2 cos(2/3) 
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(Dl z>2) basis. In this basis the squared-mass matrix reads: 

1^ ml mi + im|cos(2/3)J ^""^ 

and it may be conveniently diagonahzed by a rotation into a basis defined by the CP-even 
and CP-odd iy- sneutrino eigenstates Q: 

(r^^*;^^:! (4.4, 

The states and i'- are also mass eigenstates. The squared-mass eigenvalues are easily 
computed: 

= w-i + 2"*"^ cos(2/3) lb ml (4.5) 

which implies Am^ = m^ — m^ = 2m^. The mixing angle is fixed at the value 6 = 7r/4. We 
name this models models". They have the same field content as the minimal MSSM, 
but include the additional Jj terms. The parameters relevant for the sneutrino sector are 
mi and ms- 

The Z coupling to sneutrinos is non-diagonal in the P_) basis, which are now non- 
degenerate in mass. The first consequence is that the invisible Z-width decay occurs via 
the process Z — > viV2- By keeping mi > 100 GeV (which assures that the charged lepton 
mass bounds are satisfied) and using suitable values of the off-diagonal parameter, we 
may obtain lightest sneutrinos lighter than in the minimal MSSM models. From Eq. ( [4.5| ) 
is clear that when m_B is large (close to its allowed upper limit, which is due to preventing 
a tachionic sneutrino) mi can be small. At the same time m2 is large and the sum of the 
two never gets smaller than mz (for our parameter intervals). In this case, the Z invisible 
width does not get contributions and light sneutrinos are possible. 

In models, however, we have an additional bound, which is related to neutrino 
physics. The Jj terms may induce radiative contributions to the neutrino masses js], |l3|, 35 



28]. At 1-loop, these corrections arise from diagrams involving sneutrinos and neutralinos. 



We include this constraint by calculating the 1-loop radiative contribution to the neutrino 
mass ^fn^p!!Z\nn ^ detailed in Ref. p5|, and then imposing that |AmJ|p°°^. | does not 



"neutrino ^--^'^--•^^ ^--^ ....j^v.^.^.^, I "neutrino I 

exceed the experimental upper bound on the neutrino mass. This constraint coming from 
the 1-loop contribution ^fTL^^ntrino introduces a direct connection of the sneutrino sector 
to neutrino physics. The 1-loop correction | Am^lf^"^;^^^ | is basically proportional to the 
mass difference between the two mass eigenstates Amgneutrino = — n-i (which in turn 
depends on the parameter m^, as well as m^r). Sneutrino dark matter phenomenology will 
be therefore bounded by neutrino physics in a non trivial way. This fact will occur also for 
the class of models of Section [s], where a J term will also be present. 



In Fig. |27| we show the correlation between I ^w-ni'^^trino I ^"^^ sneutrino mass differ- 



ence, for a scan of the sneutrino mass parameters varied in the interval: 80 GeV < m/, < 
1000 GeV and 10~^ GeV < mB < 10^ GeV. In this figure, the lightest neutralino has been 
assumed to be a pure bino of 1 TeV mass. As for the neutrino mass bounds, we show in 
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Figure 27: ^ models - Absolute value of the 1-loop contribution to the neutrino mass as a 
function of the mass difference of the two sneutrinos. The sneutrino mass parameters are varied as: 
80 GeV < mi < 1000 GeV and lO"'^ GeV < ms < 10^ GeV. The lightest neutralino is a pure bino 
of 1 TeV mass. The horizontal lines denote the upper limits on the neutrino mass, as labelled. 



the plot the limits of kinematical origin: 2 eV for electron-type neutrinos |39[|, 0.2 MeV for 
muon-type neutrinos |^9| and 18 MeV for tau neutrinos [^]. Solar and atmospheric neu- 
|25|, p6|, p7| as well as cosmological bounds on massive light neutrinos 



trino results [23 
12 



104, 105] are not compatible with the kinematical bounds on muon and tau neutrinos 



listed above, unless more than 3 families are present and/or additional sterile neutrinos 
are introduced with special mixing patterns with active neutrinos. A mass bound of 2 eV 
is instead a more reliable upper bound on the neutrino mass |27, 104, 105], and our main 
conclusion will be based on this case. 



Fig. 27 shows that a mass bound of 2 eV on the neutrino mass implies a strong degen- 
eracy between the two sneutrino mass eigenstates. In this case, sneutrino phenomenology 
does not deviate significantly from the minimal MSSM discussed in Section ^. Only when 
a larger mass bound for neutrinos is allowed, the split between the two mass eigenstates 
becomes sizeable, reaching the level of 100 MeV for the muon neutrino kinematical mass 
bound, and even 30 GeV for the tau neutrino kinematical bound. These large corrections 
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Figure 28: ^ models - Absolute value of the 1-loop contribution to the neutrino mass as a function 
of the mass of lightest sneutrino. The parameters as as in Fig. |2^. The horizontal lines denote the 
upper limits on the neutrino mass, as labelled. 
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, and the ensuing large sneutrino mass splittings, occur when one of the sneu- 



trino mass eigenstate is light, opening up the possibility of light sneutrinos, which were 



precluded in the minimal MSSM. This is shown in Fig. 28, where the 1-loop correction 



I ^"^neutrino I plotted versus mi. This shows that, under this hypothesis, sneutrinos as 
light as 10 GeV are possible. 

The off-diagonal Z coupling affects also the relic abundance calculations |12|. In this 
case, the lightest sneutrino is actually co-annihilating with the heavier eigenstate, when a 
Z-boson is exchanged (while this is not the case, for instance, for higgs exchange: in that 
case the coupling is diagonal). Therefore, when the mass splitting is large the Z-mediated 
(co) annihilation cross section gets reduced and the relic abundance may increase. The relic 



abundance for the ^ models is shown in Fig. 29, for a mass bound on the neutrino mass 



of 18 MeV. This would correspond to the case of a tau sneutrino dark matter, since we 
are applying the tau neutrino kinematical mass bound. We see that the relic abundance 
increases, and light sneutrinos down to 10 GeV are acceptable. If we reduce the mass 
bound on the neutrinos down to to 2 eV, this figure reduces to the case of minimal MSSM 
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Figure 29: Jj models - Sneutrino relic abundance ^h? as a function of the sneutrino mass mi. 
The [blue] band denotes the relic abundance for allowable models due to a variation of the sneutrino 
parameters in the intervals: 80 GeV < tul < 1000 GeV and 10"'' GeV < ms < 10^ GeV. The 
other relevant supersymmetric parameters are: lightest neutralino 30% heavier than the sneutrino, 
higgs masses of 120 GeV for the lightest CP-even higgs h and 400 GeV for the heaviest CP-even 
H and the CP-odd A. The horizontal solid and dotted lines delimit the WMAP interval for cold 
dark matter. 



discussed in Section 



The off-diagonal Z coupling has impact also on direct detection |12|. In this case, the 
sneutrino-nucleus scattering for Z exchange is no longer elastic, since we have to produce 
the heavier eigenstate. When the mass difference between the two sneutrinos is large 
enough, elastic scattering in the detector is suppressed. The lightest sneutrino can scatter 
on the nucleus only when: 



Am < 



(4.6) 



2 (mi + m_v) 

This is a nice realization of inelastic dark matter, which was discussed extensively in con- 



nection with direct detection in Refs. M, 10, 11]. 



We therefore have to include this effect in the calculation of the direct detection cross 
section. This inelasticity effect produces a suppression in the direct detection rates, which 
depends on the energy of the recoil, the type of nucleus and on the energy threshold of the 
detector. We define a suppression factor for direct detection as: 



n{Ei,E2;Am) 



(4.7) 
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Figure 30: ^ models - Suppression factor of the direct detection rate due to the off-diagonal 
Z-coupling, plotted vs. the absolute value of the 1-loop contribution to the neutrino mass. The 
sneutrino mass parameters are varied as: 100 GeV < rriL < 1000 GeV and 10"'* GeV < ttlb < 
10^ GeV. [Blue] crosses refer to the Iodine nucleus, open [green] circles to the Germanium nucleus. 
The vertical lines denote the upper limits on the neutrino mass, as labelled. 



where TZ{Ei, E2; Am) denotes the direct detection rate integrated in the energy range 
{El, E2) and calculated for a sneutrino mass difference Am, while TZ{Ei,E2; 0) is the same 
quantity calculated for Am = 0, i.e. in the standard case. The quenched energies Ef Ef 
have been chosen to be: Ef = 2 KeV and — 10 KeV for the Iodine (representative for 
the DAMA/NAI experiment) and Ef' = 10 KeV and = 100 KeV for the Germanium 
(representative for the CDMS experiment). The chosen values of E^^ correspond to the 
threshold energy of DAMA/Nal and CDMS. We apply this reduction factor in the following 
way: instead of modifying the experimental result (which we cannot do separately for each 
configuration of the model parameter space), we compare the experimental results obtained 
for the standard case with a redefined scattering cross section: 



J- (scalar) 
^ '^nucleoli 



c,'t^(scalar)xZ , / ^^(scalar) x / , o\ 



eff 

where we have applied the reduction factor only to the Z-mediated cross section. These 
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Figure 31: models - Relative suppression in the direct detection rate on Germanium and 
Iodine, as a function of the sneutrino mass mi. The sneutrino mass parameters are varied as: 
100 GeV < mi < 1000 GeV and IQ-* GeV < ms < 10^ GeV. 



suppression factors (which apply to the Z-exchange case only) are shown in Fig. ^ for the 
case of a Iodine and a Germanium nucleus. The relative reduction between scattering on the 
Iodine and Germanium nuclei is shown in Fig. where it is manifest that for sneutrinos 
lighter than about 200 GeV the detection rates in CDMS can be much more suppressed than 
in the DAMA/NAI experiment. This is a consequence of the different responses of the two 
detectors to the dark matter velocity distribution function and mass, as a consequence of 
the different quenches and threshold energies. This is a practical realization of the inelastic 
dark matter scenario able to explain why current CDMS sensitivity could be insensitive to 
some cross sections which explain the DAMA/NAI effect ||9|, |l^, We see that this is 
indeed a possibility for sneutrinos in LR models. 



The effective direct detection cross section of Eq. (4.8) is shown in Fig. p3. The 



combination of the reduction effect due to the inelasticity of the Z-exchange contribution 
and the increase of the relic abundance, which affect C^'it^deon^^^^^Sh t^ie rescaling factor 
^, produces a small effect for heavy sneutrinos. On the contrary, light sneutrinos possess a 
wide range of variation of the direct detection cross section, most of which are not excluded 
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Figure 32: Jj models - Sneutrino-nucleon scattering cross section C'^nu'^loon ^ function of the 
sneutrino mass mi, for the same set of parameters of Fig. The suppression factor in the direct 
detection rate, due to the off-diagonal Z-coupling, is included in the theoretical predictions. The 
dashed, dotted and dot-dashed curves denote the CDMS upper bounds ji^, 47 1, as in Fig. ||. 



by experimental constraints. This is due to the large mass splitting of the two sneutrino 
states which effectively suppresses direct detection through Z exchange. By adopting the 
highest (dashed) line as a conservative limit from CDMS, we see that whenever mi < mz/2 
sneutrinos in ^ models are viable dark matter candidates. This nevertheless occurs if 
we allow the 1-loop corrections to the neutrino mass to be as large as the tau-neutrino 
kinematical mass limit: when a more reliable bound of 2 eV is assumed, we basically recover 
the minimal MSSM case (which also implies a lower bound on the sneutrino mass of about 
80 GeV). 



In conclusion, once a neutrino mass bound of 2 eV is assumed, the X' models do not 
exhibit a phenomenology much richer than the minimal MSSM. Direct detection almost 
excludes this possibility, except for the occurrence of mass-matching condition between the 
sneutrino mass and the Z or higgs masses (like in the case of minimal MSSM), and also in 
this case the compatibility with direct detection is marginal. 

For these reason we do not elaborate any further on the models, and we do not 
discuss neither the full super symmetric parameter-space scan neither indirect detection 
signals for this case. 
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5. Models with right— handed sneutrinos and lepton— number violating in- 
teractions: the case for a see-saw neutrino mass 



A supersymmetric model which can accommodate a Majorana mass-term for neutrinos 
and explain the observed neutrino mass pattern, and which relies on a renormalizable 
lagrangian, may be built by adding to the minimal MSSM right-handed fields and 
allowing for Jj terms . The most general form of the superpotential which accomplishes 
this conditions is [0, ^] : 

W = e^j{iiH}H] - Yi"h}l]R^ + Vj-^ HfL^jN-^) + ^M^-^N^N^ (5.1) 

where M^-^, Y/-^ and YJ'^ are matrices, which we again choose real and diagonal. For 
the Jj parameters we therefore assume: M^'^ = M 6^'^ , in order to reduce the number of 
free parameters. The Dirac mass of the neutrinos is obtained as: mj^ = V2YJ^ , while 
represent a Majorana mass-term for neutrinos. 

The general form of the soft supersymmetry-breaking potential may be written as |p5|| : 



[(ml)^-^iV^iV-^ + e^j{Kl-^ Hmk^ + Ai^Hf^N-^) + h.c] (5.2) 

where we again use the same assumptions of diagonality in flavour space for the matrices 
{Mj^y-^, {M'^y-^ , {miy'^, Aj'^ and A^-^ as we already did before. The Dirac-mass pa- 
rameter is not chosen as a free parameter: it is instead derived by the condition that the 
neutrino mass is determined by the see-saw mechanism. In this case: ml = m^y/M"^. In 
our analyses we will fix, for definiteness, the neutrino mass to be 2 eV, in order to saturate 
the bound which comes from the radiative contribution to the neutrino mass discussed in 
the previous Section. The naming convention for this class of models is "MA J models"". 

Sneutrinos now are a superpositions of two complex fields: the left-handed field vl 
and the right-handed field A^. Since we introduced Jj terms, it is convenient, as we did in 
the previous Section, to work in a basis of CP eigenstates. By defining them accordingly 
to Eq. (|4.4| ), the mass matrix for the state vector ^^j^^j = {I'X A^ P* A* ) has a the form: 

= (5.3) 
(m\ + D'^ + m\) F'^ + moM \ 

F^ + mnM m\ + M"^ + m^, + m\ 

m\ + D'^ + m\, F^ - moM 

\ F'^-mnM mj^ + M"^ + - m%J 



where = 0.5 cos(2/3) and F^ has already been defined in Eq. ( |3.4| ). The free 
parameters of the models for sneutrino sector are therefore: m^, mj\f, M, ms and F^. 

Sneutrino mass eigenstates are obtained by diagonalizing Eq. (5.4). We define them 
as follows: 

i>i = ZiiD+ + A+ + Zi3P_ + Z,4A_ i = 1, 2, 3, 4 (5.4) 
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The lightest state, which is our dark matter candidate, may now exhibit a mixing with 
the right-handed field N and the non-diagonal nature of the Z-coupling with respect of 
the CP eigenstates. These models therefore share the properties of both LR models and 
Jj models, but bring in a more rich set of parameters. In general, the terms F'^ it ttidM 
induce left-right mixing on the sneutrino eigenstates, while the m?^ term is responsible for 
CP splitting. The new Major ana-mass parameter M may lead to left-right mixing for its 
presence in the off-diagonal elements of the mass matrix A^maJ' ^ind, if it is large, can drive 
the mass of the two heavier mass eigenstates, since it enters also in the diagonal elements. 
More specifically, if = M can lead only to a left-right mixing, while if F 7^ it can 
contribute also to the CP splitting. A sizeable splitting however occurs when m£)M ~ 
are of the order of the diagonal elements. 

The actual phenomenology of the lightest sneutrino is therefore due to the relative 
values of the various parameters. While uil, which is in common with the charged slepton 
sector, has to be necessarily larger than about 80-100 GeV to fulfill the charged sleptons 
mass bounds, the other parameters are free in nature and not directly related to the 
electroweak symmetry-breaking and its scale (at least at the tree level). The right-handed 
parameter tjin and the CP-splitting parameter niB may be taken to be of the same order 
of magnitude as m^, but since they are in general independent from m/, we assume them 
to vary freely, as we did in the previous Sections. The same occurs for the parameter 
F"^. The parameter M is related to the Majorana mass of neutrinos: a natural scale for 
this parameter is therefore a high-energy scale, larger than 10^ GeV, in order to obtain a 
see-saw neutrino mass below the eV scale with a Dirac-mass term of the order of the GeV 
scale. However, a value of M = 1 TeV which produces a eV neutrino mass with mo ^ 1 
MeV is an equally viable possibility, since m^i is a Dirac-type mass which originates from 
a neutrino Yukawa interaction of unknown strength. We therefore distinguish two cases in 
our analysis: a low-scale Majorana mass, which we fix at the value of 1 TeV (this class 
of models are named for convenience MAJ[A]); a large-scale Majorana mass, which we fix 
at 10^ GeV (models MAJ[B]). We will see that a different phenomenology arises. For a 
further discussion of these parameters and their ranges, see also Ref. | p5| . 

The experimental constraints we adopt for the MAJ models (in addition to the one 
discussed in Appendix 0) are the invisible Z-width for all the relevant mass eigenstates 
at hand and the neutrino mass bound on the 1-loop corrections to the neutrino mass 
I ^ "^neutrino I — ^ ^ ■ With our choicc of the parameters, we do not have terms which mix 
different families and therefore we do not generate contributions to the flavour-violating 
decays of charged leptons, as in general may be the case |35]. We therefore do not have 
additional bounds from limits on processes like // — > e + 7 decays. 

5.1 Models with a TeV scale Majorana mass term 

Let us flrst dicuss the case of TeV-scale Majorana mass term M. In this case, seizeable 
mass splittings are possible and the mass bound from the radiative contribution to the 
neutrino mass is relevant. Fig. 33 shows I ^"^nlTtrino I versus the mass splitting of the 



neutrino I 

two lightest sneutrino eigenstate, for M = 1 TeV and a scan of the sneutrino relevant 
parameters in the ranges: 10^ GeV < mj^ < 10^ GeV, 1 GeV < uib < 10^ GeV and 
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Figure 33: MA J [A] models - Absolute value of the 1-loop contribution to the neutrino mass 
as a function of the mass difference of the two lightest sneutrino CP eigenstates, for the case of 
a Majorana-mass parameter M — 1 TeV. The other sneutrino mass parameters are varied as: 
10^ GeV < mAT < 10^ GeV, 1 GeV < tub < 10^ GeV and 1 GeV^ < < 10'' GeV". The 
horizontal lines denote the upper limits on the neutrino mass, as labelled. 



1 GeV^ < < 10^ GeV^. We see that, contrary to the case of the pure models 
discussed in the previous Section, large sneutrino mass splittings are possible even with 
a neutrino mass bound of 2 eV. This is related to the fact that the large mass splittings 
are here related to a large mixing with the right-handed field N: this therefore suppresses 
the sneutrino couplings and as a consequence also the 1-loop contribution to the neutrino 
mass. This is manifest in Fig. ^ where we show, for the same variation of parameters, 
the distribution of the sterile component of the lightest sneutrino state versus the mass 
splitting. Typically, when the sneutrino is mostly left-handed the mass splitting is also 
small, and the sterile component increases with the increase of the mass separation. 

The large sneutrino mass splitting which is now allowed also for a neutrino mass bound 
of 2 eV is relevant for the direct detection suppression effect: this is at variance with the 
]L models, where this bound implied a strong restriction of the sneutrino phenomenology, 
bringing it to the minimal MSSM case. In addition to the mass splitting, a sizeable sterile 
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Figure 34: MA J [A] models - Component of the lightest sneutrino into sterile fields as a function 
of the mass splitting of the two lightest sneutrinos, for the case of a Majorana-mass parameter 
M = 1 TeV. The other sneutrino mass parameters are varied as in Fig. |33|. 



component now also appears. The suppression factor is shown if Fig. 35, which can be 



compared to the corresponding case of fj models of Fig. 30. 

The relic abundance of the MA J [A] models is shown in Fig. It is remarkable that 
in the whole mass range from 5 GeV to 1 TeV sneutrinos can explain the required amount 
of dark matter in the Universe. For the same mass range, sneutrinos may as well be a 
subdominant component. 

Direct detection is shown in Fig. We see that three different populations arise: 
configurations on the upper right are clearly excluded by direct detection searches. Most of 
them refer to subdominant sneutrinos. Configurations on the lower right part of the plot are 
allowed but well below current direct detection sensitivity. Configurations on the center and 
left part of the plot all fall inside the current sensitivity region, and in this specific plot show 
that could explain the annual modulation effect observed by the DAMA/Nal experiment. 
We notice that a large fraction of the configurations which fall inside the DAMA/Nal 
region are also cosmologically dominant. The correlation of the direct detection cross 



section C<^^^cieon^^^ relic abundance are plotted in Fig. |38|, which shows how almost all 
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Figure 35: MA J [A] models - Suppression factor of the direct detection rate due to the off-diagonal 
Z-coupling plotted vs. the absolute value of the 1-loop contribution to the neutrino mass, for the 
case of a Majorana-mass parameter M — 1 TeV. The other sneutrino mass parameters are varied 
as in Fig. [Blue] crosses refer to the Iodine nucleus, open [green] circles to the Germanium 
nucleus. The vertical lines denote the upper limits on the neutrino mass, as labelled. 



the cosmologically relevant configurations are under investigation or under reach of direct 
detection studies. This can be confronted with the case of pure LR models shown in Fig. 



Iq , where a large fraction of the cosmologically relevant configurations were either already 



excluded by direct detection or with very low values of '^cr^u'c^eon- We remind that now 
^"^nudeon^^^*^™^ Suppression factor as discussed in Eq. (^^). 

Contrary to the case of LR models, indirect detection in not currently sensitive to the 
MAJ[A] models. Fig. 39 shows the correlation between the upgoing muon flux from the 
Earth as a function of i^^cieon' notice that all the configurations are below the current 
limits from neutrino telescopes. However, there is a large fraction of configurations which 
could be accessible by an increase of sensitivity, although the increase should be sizeable. 



Figs. 40 shows the distribution of the sterile components with respect to the sneutrino 
mass. We see that, when sneutrinos are light, they are typically right-handed, and therefore 
they do not couple to the Z boson (although they still have couplings with higgses). For 
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Figure 36: MAJ[A] models - Sneutrino relic abundance fi/i^ as a function of the sncutrino mass 
TOi , for the case of a Majorana-mass parameter AI = 1 TeV and for a full scan of the supersymmetric 
parameter space. The sneutrino parameters (other than M) are varied as follows: 10^ GeV < niN < 
10^ GeV, 1 GeV < ms < 10^ GeV and 1 GeV^ < < lO'' GeV*. All the models shown in the 
plot are acceptable from the point of view of all experimental constraints. The horizontal solid and 
dotted lines delimit the WMAP interval for cold dark matter. 



masses above 80 GeV they may exhibit both right- or left-handed behaviour. 

The distribution of the cosmologically allowed configurations in the sneutrino param- 
eter space is shown in Figs. 41 and Light states typically have small values of both 
niN and F^, and as it was already shown if Fig. |3^ they are almost all acceptable and 
currently explored by direct detection experiments. The distribution of configurations in 
these sectors of the parameter space are somewhat different from the analogous results for 
the LR models, shown in Fig. 17 and in that case, larger values of F'^ were required 
also for light sneutrinos, a difference which is traced into the new structure of the sneutrino 
mass matrix and neutrino couplings which arise here from the introduction of terms in 
the lagrangian. 

The antiproton flux at Tp = 0.23 GeV is shown in Fig. We notice that a fraction 
of the configurations with low mass are currently excluded by the BESS data. This is at 
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Figure 37: MA J [A] models ~ Sneutrino-nucleon scattering cross section CCT^u'(fi(,on ^ function of 
the sneutrino mass mi, for the case of a Majorana~mass parameter M = I TeV and for a full scan 
of the supersymmetric parameter space. Parameters are varied as in Fig. ^ [Red] crosses refer 
to models with sneutrino relic abundance in the cosmologically relevant range; [blue] open circles 
refer to cosmologically subdominant sneutrinos. The solid curve shows the DAMA/Nal region, 
compatible with the annual modulation effect observed by the experiment [^, |4^, |4^, . 

variance with the case of LR models, where instead the antiproton flux for light sneutrinos 
were all well below the BESS sensitivity. Another remarkable difference from the LR 
models is that now the configurations which are excluded by direct detection refer to 
very low antiproton fluxes, practically undetectable. These two features show that for the 
MA J [A] models antiproton searches and direct detection studies exhibit a high degree of 
complementarity, both in the ability to exclude model configurations and the prospects of 
detection. This is clearly summarized in Fig. where the high level of complementarity 
is manifest. It is also remarkable that a fraction of configurations which are currently under 
study by direct detection experiments (either just inside the CDMS sensitivity range or, 
even more interestingly, inside the DAMA/Nal annual modulation region) have a chance 
of detection by the future PAMELA and AMS flights, or by some slightely more sensitive 
and future experiments. 
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Figure 38: MAJ[A] models - Sneutrino-nucleon scattering cross section C'^'nuTioon sneutrino 
relic abundance flh^, for the case of a Majorana-mass parameter AI — 1 TeV and for a full scan 
of the supersymmetric parameter space. Parameters are varied as in Fig. The horizontal 

[green] band denotes the current sensitivity of direct detection experiments; the vertical [yellow] 
band delimits the WMAP interval for cold dark matter. 



The antideuteron searches are as well appealing, as Fig. 45 shows: the configurations 
for light sneutrinos above the GAPS sensitivity line are almost all excluded already by the 
BESS data on antiprotons, but those under the GAPS and AMS sensitivity line are under 
reach. Figs. ^ and 45 also show that for MAJ[A] models, antimatter searches are not 
sensitive to heavy sneutrinos. 

Also the gamma-ray signal is interesting, since EGRET and GLAST are sensitive to 
sneutrino dark matter also for dark matter profiles of the NEW type [106, 107 1. This is at 
variance with the LR models, where a steeper profile was needed in order to obtain sizeable 
fluxes. Fig. ^ shows the gamma-ray flux at = 1.5 GeV from the center of the Galaxy, 
for a EGRET-like angular resolution. Like for the antimatter signals, only sneutrinos with 
masses below 80 GeV produce sizeable signals. There are configuration which are at the 
level of explaining the EGRET excess. GLAST will be sensitive to a large fraction of the 
sneutino configurations. 



-52- 




]^(j-(ajQ-iHQ-w ]_Q-s ig-* jp-? jp-fi 



Figure 39: MAJ[A] models - Upgoing muon flux from sncutrino pair annihilation in the center 
of the Earth ^^^'^^^ vs. the sneutrino-nucleon scattering cross section C'^jju'^^oo'il, for the case of a 
Majorana-mass parameter M ^ I TeV and a full scan of the supersymmetric parameter space. 
Parameters are varied as in Fig. ^ [Red] crosses refer to models with sneutrino relic abundance 
in the cosmologically relevant range; [blue] dots refer to cosmologically subdominant sneutrinos. 
The horizontal solid line denotes the current upper limit from neutrino telescopes. The vertical 
solid line denotes a conservative upper limit from direct detection searches and the [green] vertical 
shaded area refers to the current sensitivity in direct detection searches. 



5.2 Models with a large scale Majorana mass parameter 

As a second example of this class of models we consider the case of a large-scale Majorana 
mass-parameter, which we fix at M = 10^ GeV. We call these models MAJ[B]. In this case, 
the large value of M drives the two highest sneutrino mass eigenstates to be very heavy 
and decoupled from the rest. The model for the light sneutrino sector is somehow less rich 
phenomenologically as compared to the case of a TeV-scale Majorana mass-parameter. 
The 1-loop neutrino mass contribution | ^"^neutrino I shown in Fig. 47 as a function of 



the sneutrino mass difference. Contrary to the case of Jj models, larger sneutrino mass 
splittings are allowed, up to almost 10 MeV. This is an interesting range for the inelasticity 
properties in direct detection dark matter [^, 10, 11]. Directly from the MA J [A] case. 
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Figure 40: MAJ[A] models - Distribution of cosmologically acceptable models in the mi vs. 
sterile-component plane, for the case of a Majorana-mass parameter M = 1 TeV and a full scan 
of the supersymmetric parameter space. Parameters are varied as in Fig. ^ [Red] crosses refer 
to models with sneutrino relic abundance in the cosmologically relevant range; [blue] dots refer to 
cosmologically subdominant sneutrinos; open dots mark the configurations which have a direct- 
detection cross section in the current sensitivity range. 



instead, larger mass splittings are not possible now. This is related to the fact that light 
sneutrinos are not obtained in this class of models. This is shown in Fig. 48 where the relic 
abundance is plotted versus mi. We see that now cosmologically relevant sneutrinosare 
present only for very restricted mass ranges, one around 80-90 GeV and the second at 
500-600 GeV. For masses above 600 GeV sneutrinos are cosmologically excluded. In the 
whole interval 90-500 GeV, sneutrinos are subdominant. 

Direct detection comes back to be a strong constraint. Fig. ^ shows that direct 
detection bounds exclude most of the configurations, and in particular all the cosmologically 
relevant ones. Nevertheless, sneutrinos in the mass range 90-300 GeV are allowed, although 
subdominant. 

Indirect detection for MAJ[B] models is not very appealing: the antiproton flux at 
Tp = 0.23 GeV is shown in Fig. |5^: configurations which pass the direct detection bound 
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Figure 41: MAJ[A] models - Distribution of cosmologically acceptable models in the rrii - 
plane, for the case of a Majorana-mass parameter M — I TeV and a full scan of the supersymmetric 
parameter space. Parameters are varied as in Fig. ^ [Red] crosses refer to models with sneutrino 
relic abundance in the cosmologically relevant range; [blue] dots refer to cosmologically subdominant 
sneutrinos; open dots mark the configurations which have a direct-detection cross section in the 
current sensitivity range. 

all refer to very low antiproton fluxes, practically undetectable. A similar situation occurs 
for antideuterons and gamma-rays. 

6. Conclusions 

In this paper we have analyzed the possibility to have sneutrinos as dark matter candi- 
dates. We have discussed various supersymmetric models, and specifically extensions of the 
Minimal Supersymmetric Standard Model which could lead to explanations of the neutrino 
mass problem. 

We have first re-analyzed the minimal version of the MSSM. Sneutrinos are here 
typically subdominant dark matter components, with low values of the relic abundance in 
all the mass range from 50-70 (their lower mass bound in the MSSM) up to 700 GeV. They 
may represent the dominant dark matter component for masses in the 600-700 GeV range. 
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Figure 42: MA J [A] models - Distribution of cosmologically acceptable models in the mi - tun 
plane, for the case of a Majorana-mass parameter M — I TeV and a full scan of the supersymmetric 
parameter space. Parameters are varied as in Fig. ^ [Red] crosses refer to models with sneutrino 
relic abundance in the cosmologically relevant range; [blue] dots refer to cosmologically subdominant 
sneutrinos; open dots mark the configurations which have a direct-detection cross section in the 
current sensitivity range. 

This possibility is actually excluded by direct detection searches, virhich allow sneutrinos to 
be a subdominant dark matter component only very marginally, and for mass- matching 
conditions that pose the sneutrino annihilation cross-section on one of the higgs poles or 
the Z pole. 

The inclusion of right-handed superfields to the supersymmetric lagrangian allows for 
a much richer phenomenology. The mixing with right-handed fields suppresses the Z~ 
coupling and leads to some increase of the relic abundance and to some decrease of the 
direct detection rate. The reduced Z coupling is also instrumental in allowing light sneutri- 
nos, by circumventing the invisible Z-width bound. From a full scan of the supersymmetric 
parameter space, we find that cosmologically dominant relic sneutrinos are present in the 
mass range from 15 GeV up to 1 TeV (where we stop our scan). When the full supersym- 
metric parameter space is considered, we find that 15 GeV is actually the mass lower bound. 
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Figure 43: MAJ[A] models - Antiproton flux at the antiproton kinetic energy Tp = 0.23 GeV as 
a function of the sneutrino mass rni, for the galactic propagation parameters which provide the 
median value of antiproton flux and for a solar activity at its minimum. The plot refers to the case 
of a Major ana-mass parameter M = 1 TeV and a full scan of the supersymmetric parameter space. 
Parameters are varied as in Fig. [Red] crosses refer to models with sneutrino relic abundance in 
the cosmologically relevant range; [blue] dots refer to cosmologically subdominant sneutrinos; light 
gray points denote configurations which are excluded by direct detection searches. The [yellow] 
shaded area denotes the amount of exotic antiprotons which can be accommodated in the BESS 
data ]f§, |§ . The dashed and dotted lines show the PAMELA [|| and AMS @ sensitivities to 
exotic antiprotons for 3 years missions, respectively. 



induced by the cosmological limit on the relic abundance. Direct detection is acceptable 
for all the allowed mass range. We also find that cosmologically dominant sneutrinos are 
accepted by direct detection, and a large fraction of the supersymmetric configurations 
predict direct detection rates at the level of the current experimental sensitivities (includ- 
ing the possibility to explain the DAMA/Nal annual modulation effect). Indirect detection 
rates offer good possibilities: antiproton fluxes are under reach of the PAMELA and AMS 
detectors in the 50-200 GeV mass range. The configurations accessible to indirect searches 
are typically cosmologically subdominant. The same occurs also for the antideuteron sig- 
nal, which is accessible by GAPS and AMS to the same configurations to which PAMELA 
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Figure 44: MAJ[A] models - Antiproton flux at the antiproton kinetic energy Tp = 0.23 GeV vs. 
the sneutrino-nucleon scattering cross section C'^nuTioon' ^^'^ case of a Majorana-mass parameter 
M = 1 TeV and a full scan of the supersymmetric parameter space. Parameters are varied as in Fig. 
|3^ . [Red] crosses refer to models with sneutrino relic abundance in the cosmologically relevant range; 
[blue] dots refer to cosmologically subdominant sneutrinos. The horizontal solid line denotes the 
upper limit from BESS |8^ and the [yellow] shaded area shows the amount of exotic antiprotons 
which can be accommodated in the BESS data. The dashed and dotted horizontal lines show the 



PAMELA |83| and AMS |84 sensitivities to exotic antiprotons for 3 years missions, respectively. 
The vertical solid line denotes a conservative upper limit from direct detection searches and the 
[green] vertical shaded area refers to the current sensitivity in direct detection searches. 



and AMS are sensitive for antiprotons. This offers a great opportunity for dark matter 
searches: a signal detectable in one antimatter channel by two different detectors, will be 
detectable also in the other channels, again by two different detectors. Gamma-rays from 
the center of the galaxy do not provide very large signals: we predict fluxes not too far 
from the EGRET excess in the 50-200 GeV mass range, but this requires a very steep 
dark matter density profile toward the galactic center (of the r~^'^ type). GLAST will be 
sensitive to configurations in this same mass range, again for a very steep profile. 

Models with a lepton-number violating term, introduced as non-renormalizable 5- 
dimensional operator, do not lead to a phenomenology very different from the standard 
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Figure 45: MA J [A] models - Antideuteron flux at the antideuteron kinetic energy (per nucleon) 
Tp — 0.23 GeV/n, as a function of the sneutrino mass mi, for the case of a Majorana-mass 
parameter M = I TeV and a full scan of the supersymmetric parameter space. Parameters are 
varied as in Fig. Notations are as in Fig. |2^. The dashed and dotted hues show the GAPS 
and AMS Q sensitivities. 



MSSM, once a mass bound of 2 eV on the 1-loop correction to the neutrino mass (which is 
induced by the terms) is considered. Only for a neutrino mass bound of 18 MeV, which 
corresponds to the kinematical mass bound for the tau neutrino, some increase of the relic 
abundance is possible. However, the direct detection limit strongly bounds these models, 
making them almost marginal. 

A renormalizable lagrangian with both right-handed fields and }L terms, which offers 
the possibility to include neutrino masses via a see-saw mechanism, again offers very rich 
sneutrino phenomenology. In the case of a TeV-scale Majorana mass-parameter, sneutrinos 
may be the dominant dark matter component for masses in the range 5 GeV up to 1 TeV 
(upper bound in our scan). The direct detection is nicely evaded in all the mass range, and 
most of the configurations fall inside the current experimental sensitivity range (including 
the possibility to explain the DAMA/Nal annual modulation effect). Antiproton fluxes 
are a stringent bound for light sneutrinos, complementary to the bound imposed by direct 
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Figure 46: MA J [A] models - Gamma-ray flux from the galactic center at the photon energy 
Ej — 1.5 GeV, as a function of the sneutrino mass mi, for a galactic profile of Moore's type |8p 
and for the angular resolution of EGRET H, ||. The plot refers to the case of a Majorana-mass 
parameter M = 1 TeV and a full scan of the supersymmetric parameter space. Parameters are varied 
as in Fig. [Red] crosses refer to models with sneutrino relic abundance in the cosmologically 
relevant range; [blue] dots refer to cosmologically subdominant sneutrinos; light gray points denote 
configurations which are excluded by direct detection searches. The [yellow] shaded area denotes 
the amount of exotic gamma-rays compatible with the EGRET excess [||, H. The dashed line 
shows the GLAST sensitivity for a 1 year data-taking and for the same EGRET angular bin. 

detection which instead is more severe for heavier particles. Many configurations for masses 
below 80-90 GeV will be explored by PAMELA and AMS, while for masses above 90 GeV 
antiproton searches loose sensitivity. Also for antideuterons, AMS and GAPS will have 
sensitivity to probe a fraction of the configurations for masses below 80-90 GeV. For 
these light sneutrinos, also gamma-rays provide a significant probe, also for NFW profiles. 
GLAST will have sensitivity to a fraction of those configurations with mass below 80- 
90 GeV. Finally, models with a large Majorana mass-parameter are strongly bounded by 
direct detection: configurations with masses in the range 90-300 GeV are not excluded 
by direct detection, but they all refer to cosmologically subdominant sneutrinos. Indirect 
detection rates are typically very suppressed. 
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Figure 47: MAJ[B] models - Absolute value of the 1-loop contribution to the neutrino mass 
as a function of the mass difference of the two lightest sneutrino CP eigenstates, for the case of 
a Majorana-mass parameter M = 10^ GeV. The other sneutrino mass parameters are varied as: 
TOAT = 0, 10^ GeV < ms < 10^ GeV and 1 GeV^ < < 10^ GeV*. The horizontal lines denote 
the upper limits on the electron neutrino mass. 

We therefore conclude that sneutrinos offer a rich phenomenology as dark matter can- 
didates, and they are a viable alternative to relic neutralinos in a wide class of supersym- 
metric models. Their phenomenology is also linked and constrained by neutrino physics 
through the problem of the origin of neutrino masses. 
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7. Appendix: The supersymmetric model 

The supersymmetric model we adopt in this paper is an effective MSSM scheme at the 
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Figure 48: MAJ[B] models - Sneutrino relic abundance flh^ as a function of the sncutrino mass 
nil, for tlie case of a Major ana-mass parameter M = 10^ GeV and for a full scan of the supersym- 
metric parameter space. The sneutrino parameters (other than M) are varied as follows: = 0, 
10^ GeV < ms < 10* GeV and 1 GeV^ < -F^ < 10'* GeV^. All the models shown in the plot are 
acceptable from the point of view of all experimental constraints. The horizontal solid and dotted 
lines delimit the WMAP interval for cold dark matter. 



electroweak scale. The free independent parameters of the model, which are not directly 
related to the sneutrino sector, are the following: the SU(2) gaugino-mass parameter M2; 
the ratio between the U(l) and SU(2) gaugino-mass parameters R = M1/M2 (in the GUT- 
induced case R = S/Stan^^^/ ~ 0.5, where 9w is the Weinberg angle); the Higgs-mixing 
parameter tan/3 = V2/V1; the mass of the pseudoscalar higgs uia', a common soft-mass 
for all the squarks mg (both right- and left-handed); a common dimensionless trilinear 
parameter for the third family A (A-^ = A^ = Arriq and Af = Amj\ the trilinear parameters 
for the other families being set equal to zero). The masses of the CP-even and charged 
higgses are calculated from ttia (and the other relevant supersymmetric parameters) by 
employing 2-loop corrections. 

The full scans of the parameter space are performed over the following ranges of 
the MSSM parameters: 1 < tan/? < 50, 100 GeV < |//| < 3000 GeV, 100 GeV < M2 < 
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Figure 49: MAJ[B] models - Sneutrino-nucleon scattering cross section ^%u'c'ioon ^ function of 
the sneutrino mass mi , for the case of a Majorana-mass parameter M = 10^ GeV and for a fuU scan 
of the supersymmetric parameter space. Parameters are varied as in Fig. ^ [Red] crosses refer 
to models with sneutrino relic abundance in the cosmologically relevant range; [blue] open circles 
refer to cosmologically subdominant sneutrinos. The solid curve shows the DAMA/Nal region, 
compatible with the annual modulation effect observed by the experiment [^, |4^, |4^, . 



3000 GeV, 100 GeV < mq < 3000 GeV, 90 GeV < < 1000 GeV, -3 < A < 3. As for 
the R parameter, we use either its mSUGRA value R = 0.5 or we scan over the interval 
0.005 < R < 0.5, depending on the case at study. In order to have the sneutrino as a dark 
matter candidate, we accept only parameter configurations for which the lightest sneutrino 
is also the lightest among all the supersymmetric particles. 

We impose the following experimental constraints: accelerators data on supersymmet- 
ric particles and Higgs-boson searches (CERN e+e~ collider LEP2 [67, 68, 54, ^] and 
Collider Detectors DO and CDF at Fermilab [ |55| , |108| ); measurements of the 6 — > s + 7 
decay process [10£]: we adopt the interval 2.89 < B(b s + 7) • 10""^ < 4.21, which is 
larger by 25% with respect to the experimental determination [109| in order to take into 
account theoretical uncertainties in the SUSY contributions [ IIC ] to the branching ratio of 
the process (for the Standard Model calculation, we employ the recent NNLO results from 
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Figure 50: MAJ[B] models - Antiproton flux at the antiproton kinetic energy Tp = 0.23 GeV as 
a function of the sneutrino mass rni, for the galactic propagation parameters which provide the 
median value of antiproton flux and for a solar activity at its minimum. The plot refers to the case 
of a Major ana-mass parameter M = 10^ GeV and a full scan of the supersymmetric parameter 
space. Parameters are varied as in Fig. [Red] crosses refer to models with sneutrino relic 

abundance in the cosmologically relevant range; [blue] dots refer to cosmologically subdominant 
sneutrinos; light gray points denote configurations which are excluded by direct detection searches. 
The [yellow] shaded area denotes the amount of exotic antiprotons which can be accommodated 
in the BESS data |7|, |§ . The dashed and dotted lines show the PAMELA l|| and AMS Q 
sensitivities to exotic antiprotons for 3 years missions, respectively. 



Ref. |111| ]); the upper bound on the branching ratio BR{Bg + n'^) |112| ]: we take 

BR{B^ fi^ + fi'^) < 1.2 • 10^'^; measurements of the muon anomalous magnetic moment 
= {g^ — 2)/2: for the deviation Aa^ of the experimental world average from the theoret- 
ical evaluation within the Standard Model we use here the range -98 < Aa^ • lO" < 565, 
derived from the latest experimental [113] and theoretical [114] data. The invisible Z-width 
constraints is also imposed on neutralinos lighter than mz/2 which occur in the gaugino 
non-universal models. 
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